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A.  SCIENTIFIC  OBJECTIVES:  The  objective  of  this  research  program  is  to  apply 
model-based  signal  processing  methods  to  enhance  the  performance  of  computation 
electromagnetics  (CEM)  simulators  for  shipboard  antenna  design.  While  recent  advances 
in  CEM  algorithms  has  significantly  reduced  the  simulation  cost  of  modeling  complex 
radiation  and  scattering  phenomena,  real-world  engineering  design  and  optimization 
often  require  that  calculations  be  carried  out  repeatedly  over  large  parameter  spaces  such 
as  frequency  and  aspect  angle,  making  the  computation  cost  still  exceedingly  high.  The 
goal  of  this  research  is  to  apply  model-based  signal  processing  algorithms  to  CEM 
simulators  to  overcome  such  computational  bottleneck  and  achieve  higher  design 
throughput.  In  particular,  we  shall  address  the  issue  of  how  to  interpolate  and  extrapolate 
the  frequency  behaviors  of  antenna  characteristics  based  on  a  sparse  set  of  computed 
data.  We  shall  also  extend  this  methodology  from  the  frequency  to  the  spatial  dimension 
to  map  out  radiation  field  behaviors  on  a  complex  platform.  Finally,  we  shall  explore  the 
transportability  of  the  methodology  to  other  applications  such  as  radar  signatures  and 
wireless  channel  characteristics. 

B.  SUMMARY  OF  RESULTS  AND  SIGNIFICANT  ACCOMPLISHMENTS:  The 

focus  during  the  first  year  of  this  research  program  is  twofold.  First,  we  set  out  to  devise 
algorithms  to  extrapolate  the  frequency  behavior  of  antenna  radiation  characteristics  on 


complex  platforms  from  a  limited  set  of  computed  data.  Our  approach  is  based  on 
applying  a  time-of-arrival  model  to  the  induced  current  on  the  platform  and  using  the 
superresolution  algorithm  ESPRIT  to  estimate  the  model  parameters.  Second,  we  set  out 
to  parameterize  complex  antenna  radiation  patterns  using  a  sparse  point  radiator  model. 
Our  approach  is  based  on  a  matching  pursuit  algorithm  that  extracts  the  relevant  radiation 
physics  contained  in  the  signal.  The  main  computational  engine  of  focus  is  the  code 
FISC,  which  employs  the  fast  multipole  technology.  The  detailed  descriptions  of  these 
two  areas  are  described  below.  They  are  followed  by  discussions  on  some  exploratory 
efforts  in  other  related  topics. 

Model-based  frequency  extrapolation  of  antenna  radiation  patterns  on  complex 
platforms.  In  antenna  design  and  analysis,  the  mounting  platform  can  have  a  significant 
effect  on  the  antenna  radiation  characteristics.  However,  rigorous  solution  of  the  radiation 
problem  over  a  complex  platform  is  very  time  consuming,  and  the  computation 
complexity  increases  dramatically  as  the  frequency  increases.  During  the  past  year,  we 
have  developed  a  model-based  frequency  extrapolation  technique  with  which  the  radiated 
field  over  a  broad  band  of  frequencies  can  be  obtained  using  rigorously  computed  results 
at  low  frequencies.  Our  approach  entails  three  steps.  First,  the  induced  current  on  the 
platform  surface  is  computed  at  low  frequencies  using  a  CEM  solver.  For  this  purpose, 
we  have  transported  the  fast  multiple  method  (FMM)  base'd  code  FISC,  originally 
developed  by  the  University  of  Illinois  on  the  UNIX  platform,  to  the  Windows  NT 
environment.  We  have  demonstrated  the  capability  of  solving  3-D  problems  of  up  to  20 
electrical  wavelengths  (over  30,000  unknowns)  on  a  desktop  PC.  We  have  incorporated 
basic  antenna  modeling  capabilities  into  FISC,  which  is  originally  developed  for  radar 
scattering  applications.  We  have  validated  FISC  results  against  measurements  and  other 
high-frequency  codes  on  several  canonical  configurations.  Lastly,  we  have  demonstrated 
the  capability  to  generate  frequency-aspect  antenna  radiation  data  for  3-D  ship-like 
platforms. 

The  extrapolation  algorithm  is  based  on  applying  the  time-of-arrival  model  to  the 
current  on  each  basis  element  on  the  surface.  In  this  model,  we  assume  that  the  total 


current  is  induced  by  different  scattering  mechanisms,  as  shown  in  Fig.  1 .  Each  of  the 
incident  mechanism  has  a  distinct  arrival  time,  so  that  the  current  can  be  written  as 

J(co)=  Zane-j(0t"  -  (1) 

n=l 

where  co  is  the  angular  frequency,  N  is  the  total  number  of  incident  mechanisms  and  t„  is 
the  arrival  time  of  the  n**1  incident  mechanism.  The  model  coefficients  a„  and  t„  are 
obtained  using  the  superresolution  algorithm  ESPRIT  [1],  which  is  based  on  the  model 
that  the  signal  consists  of  a  sum  of  exponential  and  additive  white  noise.  Given  a 
sequence  with  M  samples,  the  algorithm  can  estimate  the  number  of  exponential  N  and 
determines  the  amplitude  and  period  of  each  exponential  term.  The  basic  requirement  in 
the  number  of  samples  is  M>  2N+  1.  Once  the  model  parameters  are  found,  the  induced 
current  at  higher  frequencies  can  be  computed  using  (1).  The  radiated  field  is  then 
obtained  from  the  extrapolated  current. 


Fig.  1.  Time-of-arrival  model 


As  an  example,  we  consider  an  antenna  radiating  in  the  presence  of  a  3-D 
structure  shown  in  Fig.  2.  The  induced  current  is  rigorously  computed  over  the 
frequencies  from  0.15  to  0.24  GHz  at  an  interval  of  0.01  GHz.  The  computation  is  done 
by  FISC.  The  current  is  then  extrapolated  and  the  radiated  field  is  computed.  The 
extrapolation  result  is  compared  with  the  reference  field  data,  which  is  also  computed  by 
FISC  over  the  frequency  band  from  0.15  to  0.45  GHz  at  6  =  50°  and  0  =  0°,  as  shown  in 
Fig.  3(a).  The  corresponding  time  domain  response  obtained  via  an  inverse  Fourier 


1.5  ANTENN  A z 


Fig.  2.  Geometry  of  3-D  platform. 


(a)  Frequency  domain  response  (dB) 


Fig.  3.  Frequency  Extrapolation 


transform  of  the  frequency  data  is  plotted  in  Fig.  3  (b).  We  observe  that  the  agreement 
between  the  two  is  good  in  both  the  frequency  and  time  domains. 

While  the  above  approach  is  similar  to  that  we  had  developed  previously  in  [2] 
and  [3]  for  radar  signature  extrapolation,  several  new  challenges  have  been  identified  in 
the  present  antenna  radiation  problem.  First,  the  scattering  off  the  platform  from  antenna 
radiation  is  dominated  by  close-by  interactions  because  of  the  1/r  space  loss.  It  is 
therefore  more  difficult  to  temporally  resolve  these  closely  spaced  scattering 
mechanisms.  This  leads  to  a  large  bandwidth  requirement  in  the  CEM  computation. 
Second,  we  found  that  the  frequency  dependence  from  the  antenna  is  significant  and  must 
be  incorporated  into  the  model.  In  addition,  the  frequency  dependence  of  different 
scattering  components  of  each  current  element  could  be  different,  due  to  different 
scattering  physics.  Further  incorporation  of  these  effects  should  further  enhance  the 
accuracy  of  the  extrapolation.  Third,  we  suspect  the  results  of  FISC  for  near-field  source 
problem  may  not  be  as  accurate  as  that  from  the  far  field  incidence.  Work  is  ongoing  to 
quantify  and  further  overcome  these  limitations. 

Extraction  of  sparse  representation  of  antenna  radiation  and  coupling  data.  Task  1 
focuses  on  the  utilization  of  model-based  processing  to  extrapolate  CEM  data  for  the 
“forward”  solution  of  antenna  characteristics  in  the  presence  of  complex  platforms.  An 
equally  important  problem  from  the  design  perspective  is  the  “inverse”  algorithm  of 
spatially  pinpointing  the  locations  of  platform  scattering  based  on  CEM  data.  Toward 
this  end,  we  had  previously  developed  a  method  to  extract  sparse  model  of  the  antenna 
radiation  pattern  on  a  complex  platform  [4], [5].  This  representation  is  based  on  a  point 
radiator  model  that  describes  the  radiation  pattern  by  a  collection  of  “radiation  centers” 
on  the  platform.  The  methodology  for  obtaining  the  radiation  center  model  entails  first 
generating  the  3-D  antenna  synthetic  aperture  radar  (ASAR)  imagery  of  the  platform,  and 
then  parameterizing  the  resulting  image  by  a  collection  of  point  radiators  via  the  CLEAN 
algorithm  [6], [7],  Once  such  a  representation  is  obtained,  we  can  rapidly  reconstruct 
antenna  radiation  patterns  over  frequencies  and  aspects  with  good  fidelity.  Thus  such  a 
model  can  be  used  for  real-time  reconstruction  of  complex  antenna  patterns  in  high-level 
system  simulations.  Furthermore,  the  resulting  radiation  center  information  can  be  used 


to  pinpoint  cause-and-effect  in  platform  scattering  and  provide  important  design 
guidelines  for  antenna  placement  and  optimization.  However,  the  method  was  based  on  a 
Fourier-based  algorithm  that  relied  on  a  number  of  small-angle,  small-bandwidth 
approximations.  Furthermore,  the  concept  was  only  demonstrated  using  high-frequency 
ray-tracing  simulations  and  not  more  rigorous  CEM  data. 

During  the  past  year,  we  have  overcome  the  above  deficiencies  by  developing  a 
more  generalized  matching  pursuit  algorithm  in  the  frequency -aspect  domain  based  on  a 
point  radiator  basis.  We  have  also  demonstrated  the  feasibility  of  extracting  a  sparse 
model  of  the  antenna-platform  interaction  using  CEM  data  from  FISC.  The  matching 
pursuit  algorithm  is  implemented  based  on  the  following  radiation  center  basis: 

Es(k,d,<t>)  =  Ae-jkr0eik(xosinecos<l>+y0sin0sin<p+zocosd)  _  Ae-jkr0gj(kxx0+kyy0+kzzo) 


(2) 

where  k  is  the  wave  number,  (xo,  yo,  zo )  is  the  location  of  the  radiation  center.  The  origin 
of  the  above  basis  is  illustrated  in  Fig.  4.  Note  that  the  phase  factor  in  (2)  accounts  for 
the  path  delay  from  the  antenna  to  the  point  scatterer  on  the  platform,  and  then  to  the 
observation  direction  in  the  far  field.  To  speed  up  the  parameterization  time  of  the 
matching  pursuit  algorithm,  we  estimate  the  location  of  the  radiation  centers  by  utilizing 
the  Fourier-based  ASAR  algorithm.  The  point  with  the  highest  intensity  is  first  located  in 
the  ASAR  image  and  its  amplitude  and  coordinates  are  determined  to  serve  as  an  estimate 
of  the  strongest  radiation  center.  We  next  zoom  in  on  the  precise  location  of  the  radiation 
center  via  a  fine  search.  We  then  subtract  the  contribution  of  this  radiation  center  from 
the  total  radiated  field  and  iterate  the  process  until  the  energy  in  the  residual  signal  has 
reached  a  sufficiently  small  level. 

As  an  example,  the  radiation  center  model  is  extracted  from  the  computed 
radiation  data  from  FISC  for  the  structure  shown  in  Fig.  2.  The  frequency  range  of  the 
radiated  field  is  from  0.6  to  0.9  GHz.  The  angular  range  of  observation  is  over  a  30- 
degree  window  in  both  azimuth  and  elevation.  The  first  10  extracted  radiation  centers  are 
shown  in  Fig.  5,  with  their  amplitudes  represented  by  grayscale.  We  observe  that  the 
dominant  platform  scattering  comes  from  the  edge  diffraction  from  the  right  edge  of  the 
platform  and  the  comer  structure  formed  by  the  cylinder  and  plate.  Note  that  the 


(0,0,0) 


Fig.  4.  Radiation  center  model 


Fig.  5.  Extracted  radiation  centers  ' 


(a)  Original  radiated  field  (b)  Reconstructed  radiated  field 


Fig.  6.  Reconstruction  of  the  radiated 
field  from  the  radiation  center  model. 


radiation  center  due  to  the  edge  diffraction  is  slightly  off  the  platform.  This  is  due  to  the 
limited  resolution  of  the  matching  pursuit  algorithm.  Once  the  sparse  representation  is 
generated,  the  radiated  field  can  be  easily  reconstructed  using  the  radiation  centers.  The 
original  radiated  field  at  k2  =  0  is  plotted  in  Fig.  6(a)  as  a  function  of  the  kx  and  ky,  where 
kx,  ky  and  kz  are  defined  in  (2).  The  field  reconstructed  from  the  first  10  radiation  centers 
is  shown  in  Fig.  6(b).  The  two  patterns  match  well  over  both  frequency  and  angle.  Thus 
the  radiation  pattern  is  efficiently  approximated  by  the  radiation  center  model. 

We  shall  continue  to  refine  this  algorithm  and  examine  more  complex  topside 
structures.  In  addition,  in  preparation  for  attacking  the  near-field  antenna-to-antenna 
coupling  problem,  we  have  extended  the  ASAR  concept  to  the  near-field  coupling 
scenario  between  antennas.  This  new  algorithm  is  termed  ACSAR  (Antenna  Coupling 
Synthetic  Aperture  Radar)  [11].  We  plan  to  utilize  this  model  for  addressing  the  antenna 
coupling  issue  once  the  CEM  capability  for  simulating  antenna  coupling  characteristics  is 
in  place. 

Other  exploratory  topics.  Several  additional  studies  have  been  initiated  during  the  past 
year  and  are  described  below.  They  serve  as  possible  launching  points  into  more  relevant 
and  focused  efforts  for  the  present  program  in  the  next  year. 

We  have  begun  to  investigate  the  resonant  behaviors  of  antenna-platform 
configurations.  This  topic  is  motivated  by  the  interest  of  the  Navy  SPA  WAR  Center  in 
HF  antennas  where  ship  body  resonances  can  dominate  the  antenna  frequency 
characteristics.  In  the  lower  frequency  regime,  the  scattering  phenomenology  differs  from 
ray-optical  characteristics.  Consequently,  the  time-of-arrival  model  we  have  utilized 
successfully  may  no  longer  be  efficient  in  this  regime.  We  are  carrying  out  a  preliminary 
study  to  understand  the  scattering  phenomenology  in  these  configurations.  Once  we  have 
a  better  grasp  of  the  scattering  physics,  we  plan  to  devise  model-based  parameterization 
strategies  to  efficiently  compute  antenna  characteristics  in  the  resonant  region. 

An  important  issue  related  to  the  study  of  resonant  region  using  an  iterative  CEM 
solver  like  FISC  is  that  when  the  platform  exhibits  strong  resonance  effect,  the  iteration 
number  required  for  an  accurate  solution  can  become  large.  A  good  preconditioner  for 
the  moment  matrix  is  needed  to  alleviate  this  problem.  We  have  been  investigating  the 


use  of  wavelet  packet  basis  for  the  efficient  sparsification  of  moment  matrix  under 
previous  support  from  the  Joint  Services  Electronics  Program  [8-10].  A  possible 
extension  of  that  research  is  that  a  good  preconditioner  could  be  found  using  such  a 
wavelet  representation.  We  plan  to  further  investigate  this  topic  with  the  resonant  region 
calculation  as  a  possible  end-application. 

Finally,  we  have  begun  to  investigate  the  interpolation  problem  for  CEM  data.  In 
contrast  to  our  extrapolation  work  thus  far,  the  interpolation  problem  is  believed  to  be 
more  stable  in  general.  Although  the  potential  payoff  of  an  interpolation  algorithm  is  not 
as  great  as  extrapolation,  it  may  provide  a  more  robust  way  to  achieve  computational 
savings.  One  related  study  we  have  carried  out  recently  is  to  construct  radar  images  from 
unevenly  undersampled  data  in  the  angular  domain  [13].  We  have  devised  an  algorithm 
that  uses  a  matching  pursuit  engine  in  place  of  the  Fourier  transform  in  the  image 
construction  procedure.  The  algorithm  entails  searching  and  extracting  out  individual 
target  scattering  features  one  at  a  time  in  an  iterative  fashion.  Therefore  the  strong 
interference  between  different  target  scattering  features  can  be  reduced.  After  all  the 
main  features  are  extracted,  they  can  be  displayed  in  the  same  image  plane  without  the 
aliasing  effect  in  the  final  radar  image.  The  application  of  this  work  was  aimed  at 
constructing  images  from  measured  radar  data.  However,  based  the  success  we  have  had 
in  this  application,  we  believe  the  same  approach  may  be  applicable  for  interpolating 
CEM  data  which  are  sparsely  populated  on  an  unevenly  sampled  grid.  We  are  currently 
pursuing  this  topic  in  more  detail. 

C.  FOLLOW-UP  STATEMENT: 

During  the  coming  year,  our  research  efforts  will  be  devoted  to  two  areas.  In  the 
first  area,  we  will  continue  to  refine  the  model-based  algorithms  we  have  developed  for 
frequency  extrapolation  and  sparse  model  extraction.  The  main  objective  will  be  to  make 
the  algorithms  more  stable  and  accurate.  We  also  plan  to  study  more  complex  platforms 
in  order  to  investigate  additional  radiation  physics  and  demonstrate  the  utility  of  the 
algorithms  for  more  realistic  configurations.  In  the  second  area,  we  will  explore  several 
new  topics  outlined  in  the  last  part  of  Section  B.  We  will  begin  a  more  focused  study  of 
the  antenna  radiation  characteristics  in  the  resonant  region.  We  also  plan  to  initiate  an 


investigation  into  the  antenna  near-field  coupling  problem.  It  is  anticipated  that  some 
efforts  will  be  needed  to  upgrade  the  existing  CEM  capability  to  simulate  the  antenna 
characteristics  for  these  different  scenarios.  Finally,  we  will  extend  our  previous  work  on 
interpolating  undersampled  measured  data  to  the  present  CEM  acceleration  application. 
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ACSAR- ANTENNA  COUPLING  SYNTHETIC  APERTURE 
RADAR  IMAGING  ALGORITHM 

C.  Ozdemir  and  H.  Ling 
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Austin,  TX  78712-1084 

Abstract-A  synthetic  aperture  radar  imaging  technique  called  AC¬ 
SAR  for  antenna  coupling  scenarios  is  introduced.  It  is  shown  that  an 
ACSAR  image  of  a  platform  can  be  formed  by  inverse-Fourier  trans¬ 
forming  the  multi-frequency,  multi-spatial  coupling  data  between  two 
antennas.  Furthermore,  we  present  a  fast  ACSAR  imaging  algorithm 
that  is  specifically  tailored  to  the  shooting  and  bouncing  ray  (SBR) 
technique.  The  fast  algorithm  is  shown  to  reduce  the  total  simulation 
time  by  several  orders  of  magnitude  without  significant  loss  of  fidelity. 
Finally,  a  sparse  representation  of  ACSAR  imagery  is  introduced  by 
extracting  the  point  radiators  in  the  image.  By  parameterizing  the 
ACSAR  image,  it  is  possible  to  reconstruct  the  3-D  ACSAR  image 
and  the  3-D  frequency-spatial  data  with  a  very  sparse  set  of  radiation 
centers. 

% 

1.  INTRODUCTION 

The  electromagnetic  coupling  between  antennas  on  a  complex  platform 
is  an  important  issue  in  antenna  design.  For  an  antenna  designer,  it 
is  useful  to  identify  the  areas  on  the  platform  that  cause  strong  cou¬ 
pling  between  the  antennas.  In  this  work,  we  set  out  to  develop  an 
imaging  algorithm  to  pinpoint  the  dominant  scattering  locations  on 
the  platform  that  give  rise  to  antenna  interactions  from  the  coupling 
data  between  antennas.  Our  approach  to  this  problem  is  based  on  the 
Inverse  Synthetic  Aperture  Radar  (ISAR)  concept.  ISAR  imaging  is  a 
standard  technique  to  map  the  locations  of  dominant  scattering  off  a 
target  based  on  the  multi-frequency,  multi-aspect  backscattered  data 
[1,2] .  We  have  extended  this  concept  previously  to  the  far-field  antenna 
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Figure  1.  Comparison  of  the  ASAR  and  the  ACSAR  scenario. 


radiation  problem  by  introducing  the  Antenna  SAR  (ASAR)  imaging 
algorithm  [3].  An  ASAR  image  of  a  platform  maps  the  strong  sec¬ 
ondary  radiators  on  the  platform  from  multi-frequency,  multi-aspect 
radiation  data  in  the  far  field.  Such  information  can  be  useful  to  the 
antenna  designer  in  mitigating  platform  effects.  In  this  paper,  we  ex¬ 
tend  the  ASAR  concept  to  the  near-field  antenna  coupling  problem  by 
generating  the  Antenna  Coupling  SAR  (ACSAR)  image  of  the  plat¬ 
form  (see  Fig.  1).  By  collecting  the  multi-frequency,  multi-spatial  cou¬ 
pling  data,  it  is  shown  that  an  ACSAR  image  of  the  platform  can  be 
formed  to  display  the  dominant  scattering  locations  on  the. platform. 
To  achieve  the  required  spatial  diversity,  the  data  are  collected  on  a 

2- D  grid  at  the  receiver  site. 

This  paper  is  organized  as  follows.  First,  we  derive  the  ACSAR 
imaging  algorithm.  In  Sec.  2,  it  is  shown  that  under  the  single-bounce 
and  small-bandwidth  approximations,  a  Fourier  transform  relationship 
exists  between  the  multi-frequency,  multi-spatial  coupling  data  and  the 

3- D  locations  of  antenna-platform  interaction.  Hence,  by  3-D  inverse- 
Fourier  transforming  the  coupling  data,  it  is  possible  to  image  these 
scattering  locations  on  the  platform.  This  concept  will  be  demon¬ 
strated  by  using  the  computed  frequency-spatial  data  from  the  Shoot¬ 
ing  and  Bouncing  Ray  (SBR)  technique  [4-6].  In  section  3,  we  present  a 
fast  ACSAR  imaging  algorithm  that  is  specifically  tailored  to  the  SBR 
technique.  By  taking  advantage  of  the  ray  tracing  information  within 
the  SBR  engine,  we  demonstrate  that  a  direct  image-domain  formation 
is  possible  without  resorting  to  multi-frequency,  multi-spatial  calcula¬ 
tions.  The  image-domain  formation  algorithm  can  be  further  acceler- 
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Figure  2.  The  geometry  for  ACSAR  imaging. 


ated  by  the  application  of  the  fast  fourier  transform  (FFT)  algorithm. 
The  fast  algorithm  is  compared  to  the  general  frequency-aperture  al¬ 
gorithm  in  several  numerical  examples.  It  is  also  shown  that  the  total 
computation  time  can  be  reduced  from  hundreds  of  hours  to  minutes 
without  significant  loss  of  fidelity.  In  Sec.  4,  we  further  present  a  sparse 
model  to  represent  ACSAR  imagery.  By  parameterizing  the  ACSAR 
imagery  with  a  set  of  point  radiators  on  the  platform,  it  is  possible  to 
reconstruct  the  3-D  ACSAR  image  and  the  3-D  frequency-spatial  data 
with  a  very  sparse  set  of  radiation  centers. 

2.  ACSAR  IMAGING  ALGORITHM 

\ 

2.1  Formulation 

We  shall  first  derive  the  general  ACSAR  imaging  algorithm  that 
utilizes  multi-frequency,  multi-spatial  radiation  data  in  the  near-field 
region  of  the  platform.  We  assume  the  transmitter  and  receiver  setup 
as  depicted  in  Fig.  2.  At  the  receiver  site,  spatial  diversity  on  a  two- 
dimensional  aperture  centered  at  (a:o,0, 0)  is  used  to  achieve  resolution 
in  the  two  cross  range  dimensions.  Similarly,  to  achieve  down  range 
resolution,  frequency  diversity  is  used.  In  addition  to  the  direct  ra¬ 
diation  from  the  transmitter  to  the  receiver,  there  are  contributions 
from  the  antenna- platform  interactions.  The  scattered  electric  field  at 
the  receiver  site  due  to  the  scattering  from  a  point  P(xi,  yuZi)  on  the 
platform  can  be  written  as 
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Ea(k)  =  Ai  ■  e-jkRu  ■  e~jkR2  (1) 

where  Aj  is  the  strength  of  the  scattered  field,  Ru  is  the  path  length 
from  the  transmitter  antenna  to  P,  R2  is  the  path  length  from  P 
to  the  receiver  and  k  is  the  free-space  wave  number.  Next,  we  will 
make  two  approximations  to  the  above  equation  to  arrive  at  a  Fourier- 
based  imaging  algorithm.  The  first  assumption,  commonly  used  in 
ISAR  imaging,  is  that  the  radiation  data  are  collected  within  a  certain 
frequency  bandwidth  that  is  small  compared  to  the  center  frequency 
of  operation.  We  will  further  assume  that  the  size  of  the  aperture  at 
the  receiver  site  is  small  compared  to  path  length  R2 .  Combining 
these  assumptions,  we  can  approximate  the  phase  lag  in  the  second 
exponential  term  as  follows: 


kR2  =  kR2i  +  ko(y  cos  a*  +  2  sin  a,  sin  /?,)  (2) 

where  a*  and  fa  are  defined  in  Fig.  2,  i?2i  is  the  distance  from  P 
to  the  center  of  the  2-D  aperture,  and  the  variables  y  and  2  refer  to 
spatial  locations  in  the  receiver  aperture.  The  scattered  electric  field 
can  thus  be  approximated  by 

Es(k,  y,  z)  =  Ai  ■  •  e~jkoycosai  .  e~jk0zsina, sin/3,  ^ 

In  the  above  formula,  aTFourier  transform  relationship  exists  between 
the  variables  (k,y,z)  and  (Ri  =  Ru  +  R21,  Ui  =  ko  •  cos  ct*,  V{  = 
k0  •  sin  a{  sin  Pi) .  By  taking  the  3-D  inverse  Fourier  transform  of  the 
scattered  field  with  respect  to  k,  y,  and  z,  we  obtain  a  3-D  ACSAR 
image  of  the  platform  as  follows: 

ACSAR(R,  u,  v )  =  IFT  {Es{k,  y,  z)} 

=  IFT  {Ai  ■  e~jkR  ■  e~jko  U  ■  e~jko  V}  (4) 
=  Aj  •  S(R  -  Ri)  •  5(u  —  ui)  •  S(v  -  Vi) 

Therefore,  by  inverse  Fourier  transforming  the  multi-frequency,  multi- 
spatial  coupling  data,  the  point  scatterer  P  will  manifest  itself  as  a 
peak  in  the  image  at  (i^,  Vi)  with  amplitude  Ai .  In  practice,  since 
the  frequency  bandwidth  and  the  aperture  size  are  not  infinite,  the 
actual  spot  size  of  the  point  scatterer  will  be  inversely  proportional  to 
the  frequency  bandwidth  and  the  aperture  size.  Note  that  we  have  not 
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Figure  3.  The  geometry  for  a  multiple  bounce  mechanism. 


yet  constructed  the  ACSAR  image  in  the  original  (x,  y,  z)  coordinates. 
However,  it  is  straightforward  to  transform  the  ACSAR  image  from  the 
(i?,  u,  v)  to  the  (x,  y,  z)  coordinates  by  the  following  transformation 
formulas: 


Xq  +  R2  —  2  ■  R  ■  Xq  •  y/1  +  c 

2  •  (xo  -  R  •  Vl  +  c) 

(5a) 

(xo  -  x ) 

y  — - 1 - i— 

'  tana-  tan  p 

(56) 

z  =  tan  /?  •  (xo  —  x) 

(5c) 

o  n  1  +  tan2  Q 

c  =  tan2  (3  H - 5 — - 

tanz  a 

(5  d) 

Since  the  transformation  is  not  linear,  it  is  not  easy  to  estimate  the 
distortion  effect  of  the  transformation  on  the  ACSAR  image.  This 
issue  will  be  addressed  in  the  numerical  examples  as  well  in  Sec  4.3. 
To  summarize,  the  ACSAR  imaging  algorithm  consists  of  three  basic 
steps: 

(1)  Collect  multi-frequency,  multi-spatial  coupling  data  Es{k,y,  z) , 

(2)  Take  the  3-D  inverse  Fourier  transform  of  Es{k,y,z)  to  form 
ACSAR(i2,  u,  i/) , 

(3)  Use  (5)  to  generate  ACSAR(x,  y,  z) ,  the  image  in  the  platform 
coordinates. 

In  the  resulting  ACSAR  image,  the  dominant  scattering  points  on  the 
platform  due  to  platform  interactions  between  the  transmitter  and  the 
receiver  will  be  manifested  as  peaks. 


290 


Ozdemir  and  Ling 


Several  comments  are  in  order.  First,  the  present  Fourier-based 
imaging  algorithm  is  derived  based  on  small  bandwidth  and  small 
aperture  approximations.  Since  these  assumptions  are  similar  to  those 
made  in  the  ISAR  imaging  algorithm,  no  attempts  have  been  made  by 
us  in  this  paper  to  quantify  the  errors  related  to  these  approximations. 
Second,  the  present  ACSAR  algorithm  is  based  on  the  single-bounce  as¬ 
sumption.  Therefore,  multiple-bounce  mechanisms  will  not  be  mapped 
correctly  in  the  resulting  ACSAR  image.  This  situation  is  similar  to 
the  standard  ISAR  algorithm.  We  shall  discuss  how  the  multi-bounce 
mechanisms  are  manifested  in  the  ACSAR  image  to  guide  us  in  im¬ 
age  understanding.  We  consider  an  n -bounce  scattering  mechanism 
as  shown  in  Fig.  3.  The  corresponding  scattered  electric  field  for  this 
mechanism  is  given  by 

Es(k,  y,z)  =  A  -  e-Jfc(triptotai+#2n)  .  Q-jko-y  cosan  .  g-jko-zsina-n  sin  0n  ^ 

where  triptotal  =  trip(  + . . .  +  tripn_1  corresponds  to  the  path  that  the 
wave  travels  from  the  transmitter  to  the  last  hit  point  and  i?2n  is  the 
path  length  from  the  last  hit  point  to  the  center  of  the  aperture  at  the 
receiver  site.  By  taking  the  3-D  inverse-Fourier  transform  of  the  above 
equation,  we  obtain  a  peak  in  the  ACSAR  image  at  ( Ry, ,  un,  vn)  where 
Rn  =  triptotal  +  i?2n>  un  =  ko  •  cos  on  and  vn  =  ko  •  sino;n  •  sin/?„  . 
We  observe  that  this  corresponds  nearly  to  the  last  hit  point  on  the 
platform.  The  u  and  v  values  (and  the  corresponding  a  and  (3  an¬ 
gles)  correspond  exactly  to  those  of  the  last  hit  points,  while  the  R 
value  corresponds  to  the  cumulative  delay  of  the  mechanism.  This 
means  that  a  multiple-bounce  mechanism  will  be  mapped  in  the  AC-* 
SAR  image  as  a  point  shifted  from  the  last  hit  point  on  the  platform 
in  a  direction  opposite  to  i?2n  >  the  vector  from  the  last  hit  point  to 
the  receiver.  This  effect  will  be  seen  in  the  example  next.  Finally,  in 
the  ACSAR  algorithm  derivation  only  the  scattered  field  is  considered. 
If  the  total  field  is  applied  to  the  algorithm,  the  direct  radiation  from 
the  transmitter  to  the  receiver  will  be  imaged  at  the  transmitter  loca¬ 
tion.  Since  this  peak  will  likely  be  the  strongest  one  in  the  image,  it 
may  pose  a  dynamic  range  problem  as  the  other  scattering  mechanisms 
due  to  antenna-platform  interactions  may  be  overshadowed.  For  sim¬ 
ulation  data,  the  isolation  of  the  scattered  field  contribution  from  the 
direct  antenna  radiation  is  straightforward.  For  measurement  data,  a 
separate  calibration  run  to  measure  the  direct  antenna  radiation  may 
be  required  to  remove  its  contribution  from  the  total  field. 
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2.2  Numerical  Example 

The  algorithm  presented  above  is  first  tested  using  a  simple  config¬ 
uration  shown  in  Fig.  4.  The  platform  is  a  perfectly  conducting  square 
plate  of  size  10m  by  10m.  The  transmitter  antenna  is  a  half-wave 
dipole  (at  10  GHz)  and  is  located  20m  above  the  plate.  A  total  of  32 
frequencies  are  computed  within  a  0.35  GHz  bandwidth  at  the  center 
frequency  of  10  GHz.  The  simulation  is  carried  out  by  a  modified  ver¬ 
sion  of  the  SBR  code  Apatch  [6].  The  scattered  field  is  collected  on  an 
aperture  centered  at  (20m,  0,  0).  The  collection  aperture  consists  of 
16  x  32  =  512  points,  which  range  from  -0.66m  to  0.56m  in  16  steps 
along  the  y  direction  and  from  -0.46m  to  0.44m  in  32  steps  along  the 
z  direction.  In  the  computed  data  only  the  scattered  field  from  the 
platform  is  considered  and  the  primary  radiation  due  to  the  antenna  is 
not  included.  Using  the  ACSAR  algorithm,  we  first  generate  the  3-D 
ACSAR  images  in  ( R ,  u,  v )  domain.  While  applying  the  algorithm,  a 
3-D  Hanning  window  is  used  prior  to  the  FFT  operation  to  suppress 
the  sidelobes  in  the  image.  For  display  purpose,  the  3-D  ACSAR  im¬ 
age  is  projected  onto  the  2-D  R-u ,  R-v  and  u-v  plane  as  shown  in 
Fig.  5(a).  To  view  the  ACSAR  image  in  the  platform  (x,y,z)  coordi¬ 
nates,  we  perform  the  necessary  transformations  in  (5).  The  resulting 
3-D  ACSAR  image  is  then  projected  onto  the  2-D  y-z ,  x-y,  and  x-z 
planes  as  shown  in  Fig.  5(b).  The  dynamic  range  of  the  images  is  cho¬ 
sen  to  be  35  dB.  According  to  the  geometry,  we  expect  a  peak  in  the 
(x,  y,  z)  image  at  (10m,  —20m ,  0)  corresponding  to  the  specular  point 
on  the  plate.  The  corresponding  R,  u  and  v  values  for  this  point  are 
(44.72m,  186.4,  0).  We  observe  from  Figs.  5(a)  and  5(b)  that  the  peak 
occurs  at  the  expected  location.  However,  we  notice  that  the  specular 
peak  is  not  as  well  focused  in  the  xyz- plane  as  the  Ruv-plane .  This 
degradation  is  due  to  the  additional  (R,  u,  v)-to-(x,  y,  z)  transforma¬ 
tion. 

As  the  second  example,  we  used  a  much  more  complex  test  platform 
called  ‘Slicy’,  whose  CAD  geometry  is  shown  in  Fig.  6.  Slicy  contains 
a  number  shapes  on  top  of  the  platform  including  a  closed  cylinder, 
an  open  cylinder,  a  corner  reflector  and  a  step  region.  The  platform 
is  assumed  to  be  perfectly  conducting.  A  half-wave  dipole  is  used  as 
the  transmitter  and  is  placed  at  the  origin.  A  total  of  32  frequencies 
are  computed  within  a  0.19  GHz  bandwidth  at  the  center  frequency 
of  10  GHz.  The  scattered  field  is  collected  on  an  aperture  centered 
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Figure  6.  The  geometry  of  the  test  example  Slicy. 


at  (35m,  0, 0) .  There  are  32  x  32  =  1024  spatial  points  on  the  aper¬ 
ture,  ranging  from  —1.14  m  to  1.14  m  along  the  y  direction  and  from 
—0.49  m  to  0.49  m  along  the  z  direction.  The  total  computation  time 
for  the  SBR  simulation  is  about  160  hours  on  an  IBM  RS6000/590 
workstation.  Using  the  ACSAR  algorithm,  we  generate  the  3-D  AC- 
SAR  image  for  the  configuration.  Again,  a  3-D  Hanning  window  is  used 
to  suppress  sidelobes.  The  3-D  ACSAR  image  projected  onto  the  2-D 
R-u,  R-v  and  u-v  planes  are  shown  in  Fig.  7(a).  Fig.  7(b)  shows  the 
projected  ACSAR  image  in  the  platform  ( x ,  y,  z)  coordinates  along  the 
y-z ,  x-y  and  x-z  planes.  Also  plotted  in  the  images  is  the  platform 
overlay.  The  dynamic  range  of  the  images  is  35dB.  We  observe  that 
the  dominant  scattering  is  from  the  middle  of  the  platform  around  the 
point  (17.5m,  —40m,  0).  This  is  the  expected  specular  point  on  the 
platform  (labeled  as  mechanism  (i)).  Also  apparent  is  the  scattering 
off  the  curved  region  between  two  steps  (mechanism  (ii)).  This  mecha¬ 
nism  appears  to  be  quite  diffused  in  the  image.  We  believe  this  is  due 
to  the  distortion  effect  from  the  transformation  in  (5).  A  technique 
to  alleviate  this  effect  will  be  discussed  in  Sec.  4.3.  We  also  observe 
some  scattering  mechanisms  from  the  open  cylinder  (mechanism  (iii)). 
Since  the  open  cylinder  is  expected  to  create  multiple  bounces,  the 
contribution  in  the  image  is  shifted  in  the  direction  away  from  the  re- 
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(iii)  (i) 


Figure  7.  2-D  projected  ACSAR  images  for  Slicy 

(a)  In  v)  coordinates. 

(b)  In  (x,  y,  z)  coordinates. 


ceiver  as  discussed  in  the  last  section.  As  a  result,  it  is  located  around 
the  outside  region  of  the  open  cylinder.  Finally,  the  scattering  off  the 
top  of  the  closed  cylinder  is  labeled  as  mechanism  (iv),  although  it  is 
rather  weak. 

3.  FAST  ACSAR  IMAGE  GENERATION  USING  SBR 

The  ACSAR  image  formation  algorithm  discussed  in  the  last  section 
can  be  applied  to  either  measurement  or  simulation  data.  In  this  sec¬ 
tion,  we  shall  derive  a  fast  ACSAR  imaging  algorithm  specially  tailored 
to  the  SBR  technique.  This  is  accomplished  by  utilizing  the  ray  in¬ 
formation  within  the  SBR  ray-tracing  engine  to  generate  the  ACSAR 
image  without  resorting  to  any  multi-frequency,  multi-spatial  calcula¬ 
tions.  We  first  derive  an  image-domain  formula  of  ACSAR  imaging  [7]. 
Next,  we  implement  a  fast  ray  summation  scheme  [8,9]  based  on  the 
image-domain  formulation  to  reduce  the  total  computation  time.  It  is 
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demonstrated  that  an  ACSAR  image  can  be  simulated  at  a  fraction  of 
the  time  as  the  brute-force  frequency-aperture  approach. 

3.1  Image-Domain  Formulation  of  ACSAR  Imaging 

In  applying  the  SBR  technique  to  the  antenna-platform  interaction 
problem,  rays  are  first  shot  from  the  phase  center  of  the  transmitter 
and  traced  according  to  geometrical  optics.  At  the  exit  point  of  each 
ray  before  they  leave  the  platform  altogether,  a  ray-tube  integration 
is  carried  out  to  find  the  contribution  of  each  ray  to  the  field  at  the 
receiver.  Under  this  construct,  the  scattered  near  field  can  be  written 
as: 

Es(k,  ^  •  Ci  •  (AAU%  9)  •  e~-?k  rA  (7) 

i  rays 

where  rA  is  the  position  vector  from  the  last  hit  point  A  to  the 
receiver,  (A,4)exit  is  the  cross  section  of  the  exit  ray  tube,  S{<f>,9) 
is  the  normalized  radiation  pattern  from  the  ray  tube  and  Ci  is  the 
field  at  the  exit  ray  tube.  For  sufficiently  small  ray  tubes,  S(<j>,  6)  can 
be  assumed  to  be  unity.  In  the  image-domain  formulation,  we  will 
assume  that  the  frequency  bandwidth  is  small  compared  to  the  center 
frequency  and  the  aperture  size  at  the  receiver  site  is  small  compared 
to  the  path  length  from  the  last  hit  point  to  the  receiver.  With  these 
assumptions,  we  have  k  •  rA  =  k  •  i?2  +  kocosa  •  y  +  ko  ■  sin » sin /3  •  z . 
We  further  assume  that  the  platform  is  perfectly  conducting.  (For  non- 
perfectly  conducting  platforms,  this  derivation  still  applies  as  long  as 
the  frequency  bandwidth  is  small.)  Using  these  assumptions,  we  can 
let 

^■Ci-  (A.4)exit  »  (T{  •  e-iMtrip)*  (8) 

where  cq  is  assumed  to  be  independent  of  frequency  and  (trip),  is  the 
total  path  that  the  i'  th  ray  traveled  from  the  transmitter  to  the  last 
hit  point.  Consequently,  the  scattered  electric  field  at  the  receiver  can 
be  written  as 

Es(k,y,z)  =  •  e~jk{tlip'+R2i)  •  e-^o-cosai  y  .  e-jka  smaxsmPi  Z^ 

i  rays 

or  by  the  change  of  variables  similar  to  that  in  Sec.  2,  we  have  simply 
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Es(k,  y,  z)  —  e~jkRi  •  e~jUi  V  •  e~jV'  Z  (9*0 

i  rays 

Substituting  the  above  equation  to  our  ACSAR  imaging  formula,  we 
get 

ACSAR(R,  u,  v )  =  IFT3{£'s(fc,  y,  z)} 

=  <Ti  ■  IFT{e~jkR'}  ■  IFT{e~jvUi}  ■  IFT{e~jzVi} 

i  rays 

(10) 

The  inverse-Fourier  transform  operations  in  (10)  can  be  performed  in 
closed  form  with  the  result 

ACSAR(R,  u,  v)  =  Y;  '  i2Ak  ■  ejko<-R~Ri)  sinc(A k(R  -  Ri ))} 

i  rays 

•  {2ko&y  ■  sinc(Ay(u  —  Uj))} 

•  {2koAz  •  sinc(A2(v  —  Uj))} 

(11) 

where  A  k  is  the  half  bandwidth  in  frequency,  A y  and  A z  are  the 
half-lengths  of  the  aperture  in  y  and  z ,  respectively.  We  observe  that 
the  image-domain  formulation  gives  the  contribution  of  each  ray  to  the 
total  ACSAR  image  explicitly.  Since  cr* ,  ( Xi,yuZ{ )  are  available  from 
the  ray  tracing,  it  is  straightforward  to  calculate  the  corresponding 
(Ri,Ui,Vi)  values.  Therefore,  given  (<Tj,®t,  y^Zi)  for  each  ray,  we  can 
form  the  ACSAR  image  directly  by  summing  up  the  weighted  3-D  sine 
functions  in  (11). 

\ 

3.2  Fast  Algorithm 

The  image-domain  expression  of  (11)  is  time  consuming  to  carry 
out.  We  next  apply  a  fast  ray  summation  scheme  to  speed  up  the 
calculation  in  (11). We  rewrite  the  image  domain  formula  as 

ACSAR(i?,  u,  v )  =  ^  pi  ■  h(R  —  Ri ,  u  —  ui,v  —  Vi )  (12a) 

i  rays 

where  pi  =  8  •  cr*  •  k$Ak  ■  Ay  •  A z  and  h(R,  u,  v )  is  the  3-D  ray-spread 
function  given  by 

h(R,  u,  v)  =  eikoR  sin  c(AA:  •  R)  •  sin  c(A y  •  u )  •  sin  c(Az  •  v)  (126) 
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We  can  cast  (12a)  as  a  convolution  between  a  3-D  weighted  impulse 
function  and  h  as  shown  below: 


ACSAR(R,  u,  v) 


Pi  •  6{R  -  Ri,u-Ui,v  -  Vi)  > 
i  rays  J 


*  h(R,  u,  v ) 


To  speed  up  the  computation  time  of  the  convolution,  we  use  an  FFT- 
based  fast  scheme  proposed  by  Sullivan  [9].  Since  the  3-D  weighted 
impulse  train  g  is  not  uniformly  spaced,  we  first  interpolate  the  func¬ 
tion  onto  a  uniformly  sampled  grid.  To  ensure  the  accuracy  of  com¬ 
putation,  we  oversample  the  grid  at  N  times  the  Nyquist  rate.  For 
the  interpolation,  we  use  a  linear  approximation  to  update  the  closest 
eight  neighbors  depending  on  their  distances  from  the  location  of  the 
original  impulse.  The  detail  of  this  scheme  can  be  found  in  [3]  and 
will  not  be  repeated  here.  After  applying  the  interpolation  scheme, 
the  ACSAR  image  can  be  generated  by 


ACSAR(R,  u,  v )  =  IFFT3(FFT3(5S(R,  u, u))  • FFTz(h(R ,  u, v))}  (14) 

where  FFT3  and  IFFT3  are  the  three-dimensional  forward  and  in¬ 
verse  fast  Fourier  transforms,  respectively.  The  final  ACSAR(x,  y,  z ) 
image  can  be  generated  using  the  transformations  in  (5). 


3.3  Numerical  Example 


The  fast  approach  presented  above  is  applied  to  the  same  Slicy  ge¬ 
ometry  that  was  used  in  the  frequency-aperture  approach.  In  the  in¬ 
terpolation  grid,  an  over  sampling  ratio  of  N  =  4  is  used  in  all  three 
dimensions  to  ensure  the  accuracy  of  the  computation.  The  results  are 
depicted  in  Figs.  8(a)  and  8(b)  as  2-D  projected  ACSAR  images  in  the 
three  principal  planes.  Fig.  8(a)  shows  the  images  in  the  ( R ,  u,  v)  coor¬ 
dinates,  while  Fig.  8(b)  shows  the  images  in  the  ( x ,  y,  z)  coordinates. 
Again,  a  3-D  Hanning  window  is  applied  and  the  dynamic  range  of 
the  images  is  chosen  to  be  35dB.  By  comparing  these  images  with  the 
images  that  are  generated  by  using  the  frequency-aperture  approach 
(Figs.  7(a)  and  7(b)),  we  observe  a  fairly  good  agreement  where  all 
the  key  features  are  correctly  simulated.  The  rms  error  between  the 
frequency-aspect  result  and  the  image-domain  result  in  the  (R,u,v) 
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Figure  8.  2-D  projected  ACSAR  images  for  Slicy  using  the  fast  algo¬ 
rithm. 

(a)  Image  formed  in  (i?,  u ,  v)  coordinates. 

(b)  Image  transformed  to  (a,  y,  z )  coordinates. 

(c)  Image  formed  directly  in  (x,y,z)  coordinates. 


domain  is  calculated  to  be  2.13%.  The  image-domain  results  lead  to  a 
more  focused  image  due  to  its  small  bandwidth  approximation.  This  is 
consistent  with  our  past  experience  from  ISAR  and  ASAR  algorithms. 
To  show  the  time  savings  that  we  gain  by  using  the  image-domain 
approach,  the  total  computation  times  of  the  frequency- aperture  algo¬ 
rithm  and  the  fast  algorithm  are  compared  in  Table  1.  The  timings 
are  divided  into  the  ray-tracing  time  and  the  ray-summation  time.  As 
it  can  be  seen  from  the  table,  the  ray-tracing  time  of  9  minutes  is  the 
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Table  1.  ACSAR  simulation  time  for  Slicy  on  an  IBM  RS6000/590 
workstation. 


Ray  Tracing 
Time 

Ray  Summation 
Time 

Total 

Time 

Frequency- 

Aperture 

Approach 

'9  minutes 

161  hours 

161  hours 

Fast  Algorithm 

9  minutes 

6  minutes 

15  minutes 

same  for  both  approaches.  Ray  summation  takes  over  160  hours  for 
the  frequency-aperture  approach  and  is  the  dominant  portion  of  the 
total  simulation  time.  The  fast  approach  reduces  the  ray  summation 
time  to  only  6  minutes. 

In  addition  to  large  savings  in  simulation  time,  it  is  also  possible  to 
circumvent  the  image  distortion  effect  due  to  the  coordinate  transfor¬ 
mation  by  using  the  image-domain  formulation.  This  is  accomplished 
by  generating  the  ACSAR  image  directly  in  the  platform  (x,  y,  z)  co¬ 
ordinates  in  (12)  instead  of  the  (i2,  u,  v)  coordinates.  To  find  the 
correct  {x^y^Zi)  position  for  a  multi-bounce  ray,  we  move,  along  the 
direction  opposite  to  the  vector  from  the  last  hit  point  to  the  receiver. 
The  amount  of  the  shift  is  given  by  the  cumulative  trip  of  the  ray. 
This  has  been  discussed  earlier  in  the  image  interpretation  context 
in  Sec.  2.1.  In  this  manner,  the  image  in  the  platform  coordinate  is 
updated  one  ray  at  a  time  using  the  3-D  ray-spread  function  in  the 
(x,  y,  z)  coordinate  system,  thereby  eliminating  the  need  to  perform 
the  transformation  in  (5).  Fig.  8(c)  shows  the  result  of  generating  the 
image  directly  in  the  (x,y,z)  coordinates.  We  observe  that  all  the 
features  are  now  much  more  focused  than  those  in  Fig.  8(b). 

To  summarize,  we  have  demonstrated  the  feasibility  to  simulate  AC¬ 
SAR  imagery  using  the  SBR  technique  very  rapidly.  This  is  achieved 
by  a  combination  of  the  image-domain  formulation  and  an  FFT-based 
fast  algorithm.  Furthermore,  we  have  shown  that  it  is  possible  to  form 
the  image  in  the  platform  coordinates  directly  by  using  this  approach. 
This  eliminates  the  need  to  perform  the  required  coordinate  transfor¬ 
mation  which  can  defocus  the  features  in  the  image. 
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4.  A  SPARSE  REPRESENTATION  OF  ACSAR  IMAGERY 
4.1  Point  Radiator  Model  of  ACSAR  Imagery 

In  simulating  the  ACSAR  image  of  the  platform  using  SBR,  tens  of 
thousands  of  rays  are  shot  onto  the  platform.  Yet,  as  we  can  see  from 
the  resulting  data,  the  ACSAR  images  are  in  fact  quite  sparse.  This  is 
because  the  rays  interfere  with  each  other  to  give  rise  to  strong  coherent 
scattering  over  a  small,  localized  area  on  the  platform.  Hence,  it  may 
be  possible  to  accurately  represent  the  ACSAR  image  by  a  finite  set 
of  point  “radiation  centers”  on  the  platform.  This  concept  is  similar 
to  the  scattering  center  representation  in  radar  scattering  problems 
[10,11].  We  shall  adopt  a  point  radiator  model  which  has  the  same 
form  of  (12a)  to  parameterize  the  ACSAR  image 

ACSAR(x,  y,z)=  ^  An  ■  h(x  -  xn,y  -  yn,  z  -  zn)  (15) 

n  radiation 
centers 


where  h(x,  y,  z)  is  given  in  (12b)  and  An  is  the  strength  of  the  point 
radiator  located  at  (xn,yn,zn) .  The  remaining  problem  is  to  deter¬ 
mine  (An,xn,yn,  zn)  to  find  a  radiation  representation  for  the  image. 

4.2  Extraction  Using~CLEAN 

To  perform  the  extraction  process  in  (15),  we  apply  the  image^ 
processing  algorithm  CLEAN.  CLEAN  is  a  well-known  extraction  algo¬ 
rithm  first  developed  in  radio  astronomy  [12],  and  has  been  successfully 
applied  for  scattering  center  extraction  [10].  It  is  a  robust  technique 
that  successively  picks  out  the  highest  point  in  the  image  and  removes 
its  point-spread  function  (assumed  to  be  the  ray-spread  function  in  our 
case)  from  the  image.  At  the  n  th  iteration  of  CLEAN,  if  the  strength 
of  the  point  scatterer  is  An  ,  the  3-D  residual  image  can  be  written  as 

[3  —  D  Residual  Image]n+i  =[3  -  D  Residual  Image]n 

-  [An- h{x  -  xn,y -yn,z  -  zn)}  K 

The  above  iteration  process  continues  until  the  highest  point  in  the 
residual  image  falls  below  to  a  user-defined  threshold.  Once  we  extract 
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Figure  9.  (a)  The  locations  of  the  150  extracted  radiation  centers, 
(b)  2-D  projected  ACSAR  images  reconstructed  by  using  the  extracted 
radiation  centers. 


the  radiation  centers  from  the  image,  we  can  reconstruct  the  multi- 
frequency  and  multi-spatial  data  by  the  formula: 

Es(k,y,z)=  An  •  e~jk  Rn  ■  e~jUn  V  -  e~jv”  z  (IT) 

n  radiation 
centers 

where  ( Rn,un,vn )  axe  related  to  (xn,yn,zn)  through  the  transforma¬ 
tion  described  earlier. 

We  apply  the  CLEAN  extraction  algorithm  to  the  ‘Slicy’  example. 
We  extract  a  total  of  150  point  radiators  from  the  3-D  ACSAR  image 
shown  in  Fig.  8(c).  The  locations  of  the  extracted  radiation  centers  are 
projected  on  three  different  2-D  planes  and  are  shown  as  small  circles 
in  Fig.  9(a).  As  it  can  be  seen  from  the  figures,  most  of  the  radiation 
centers  are  located  on  the  platform  at  y  =  —40m .  The  remainders  are 
located  on  the  top  of  tall  cylinder  and  around  the  step  region.  We  then 
reconstruct  the  ACSAR  image  using  the  radiation  centers  and  they  are 
shown  in  Fig.  9(b).  We  use  the  same  dynamic  range  of  35dB  when 
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-  Original 

- Reconstructed 

Figure  10.  Comparison  of  the  original  and  the  reconstructed  field 
patterns. 

(a)  Frenquency  sweep  at  (35,  0,  0). 

(b)  Spatial  comparison  at  10  GHz  along  (35,  y,  0). 

(c)  Spatial  comparison  at  10  GHz  along  (35,  0,  z).  * 


plotting  the  images.  By  comparing  the  reconstructed  ACSAR  images 
to  the  reference  image  in  Fig.  8(c),  we  see  that  the  reconstruction 
is  good.  The  rms  error  between  the  original  and  the  reconstructed 
ACSAR  images  is  calculated  to  be  less  than  0.1%.  Next,  we  reconstruct 
the  frequency  and  the  spatial  data  by  using  the  formula  in  (17).  The 
original  (solid)  and  the  reconstructed  (dashed)  patterns  are  shown  in 
Fig.  10.  In  Fig.  10(a),  the  frequency  sweep  comparison  is  made  at 
the  center  point  of  the  receiver.  In  Fig.  10(b)  and  10(c),  the  spatial 
variation  comparisons  are  made  at  the  10  GHz  center  frequency  along 
the  y-axis  and  along  the  z-axis ,  respectively.  We  observe  a  fairly 
good  agreement  for  the  frequency-spatial  reconstruction  by  using  the 
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radiation  center  model.  Therefore,  we  can  represent  complex  coupling 
data  between  antennas  on  a  platform  by  using  only  a  small  set  of 
point  radiators.  Such  sparse  data  representation  may  be  useful  for 
data  compression  and  system  simulation  applications. 

4.3.  Eliminating  Coordinate  Transform  Distortion  Using 
CLEAN 

We  now  describe  a  method  that  takes  advantage  of  the  CLEAN  al¬ 
gorithm  to  eliminate  the  image  distortion  resulting  from  the  ( R ,  u ,  v)  - 
to-  (a;,  y,  z)  coordinate  transformation  in  the  frequency-aperture  image 
formation  algorithm  described  in  Sec.  2.  In  Sec.  3.3,  we  have  already 
shown  that  this  can  be  accomplished  by  using  the  image-domain  for¬ 
mulation.  However,  the  approach  is  only  applicable  when  coupled  to 
the  SBR  simulation.  When  data  from  measurement  or  other  types  of 
computational  electromagnetics  simulators  are  used,  we  must  adopt  a 
different  strategy.  Our  approach  is  to  first  parameterize  the  image  in 
the  (/?,  u,  v)  coordinate  system  using  the  CLEAN  algorithm.  Once  the 
image  is  parameterized  by  radiation  centers  located  at  (Rn,  un,  vn) ,  we 
use  (5)  to  find  the  corresponding  (xn,yn,2n)  location  of  each  radia¬ 
tion  center  in  the  platform  coordinate.  Then,  the  construction  of  the 
image  in  the  platform  coordinate  is  carried  out  using  (15).  Note  that 
the  main  advantage  of  this  approach  is  that  the  point-spread  response 
in  the  final  platform  image  is  well-controlled.  Consequently,  the  dis¬ 
torted  point-spread  response  in  the  (x,y,  z)  coordinates  is  avoided  in 
the  platform  image. 

As  an  example,  we  apply  this  methodology  to  the  ACSAR  image  in 
Fig.  7(a)  by  extracting  a  total  of  150  radiation  centers.  After  trans¬ 
forming  these  radiation  centers  into  the  platform  coordinates,  we  form 
the  ACSAR  image  directly  in  the  (x,y,z)  coordinates.  The  result 
is  shown  in  Fig.  11.  By  comparing  the  new  images  to  the  images  in 
Fig.  7(b),  we  observe  that  the  mechanisms  are  much  better  focused 
using  this  approach.  It  is  therefore  the  preferred  step  to  generate  the 
final  platform  image  using  the  general  frequency-aperture  algorithm. 


5.  CONCLUSION 

The  ACSAR  (antenna  coupling  synthetic  aperture  radar)  imaging  al¬ 
gorithm  has  been  introduced  to  map  platform  scattering  locations  in 
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Figure  11.  2-D  projected  ACSAR  images  based  on  frequency-aperture 
data,  formed  by  first  parameterizing  the  ( R,u,v )  image  using  CLEAN 
and  reconstructing  the  image  in  the  (x,  y,  z )  platform  coordinates. 


antenna  coupling  scenarios.  It  is  shown  that  under  the  single-bounce 
and  small-bandwidth  approximations,  a  Fourier  transform  relationship 
exists  between  the  multi-frequency,  multi-spatial  coupling  data  and  the 
3-D  locations  of  antenna-platform  interaction.  Hence,  by  3-D  inverse- 
Fourier  transforming  the  coupling  data,  it  is  possible  to  image  the  scat¬ 
tering  locations  on  the  platform.  This  concept  has  been  demonstrated 
by  using  the  computed  frequency-spatial  data  from  the  SBR  technique. 
However,  the  algorithm  is  general  and  can  be  applied  to  other  simu¬ 
lation  data  or  even  measurements.  Furthermore,  we  have  presented  a 
fast  ACSAR  imaging  algorithm  that  is  specifically  tailored  to  the  SBR 
technique.  By  taking  advantage  of  the  ray  tracing  information  within 
the  SBR  engine,  we  have  demonstrated  that  a  direct  image-domain  for¬ 
mation  is  possible  without  resorting  to  multi-frequency,  multi-spatial 
calculations.  The  image-domain  formation  algorithm  can  be  further  ac¬ 
celerated  by  the  application  of  the  FFT  algorithm.  The  fast  algorithm 
has  been  shown  to  reduce  the  total  simulation  time  from  hundreds  of 
hours  to  minutes  without  significant  loss  of  fidelity.  Finally,  we  have 
presented  a  method  to  extract  a  sparse  point-radiator  model  to  rep¬ 
resent  ACSAR  imagery.  This  is  accomplished  by  parameterizing  the 
ACSAR  image  with  a  set  of  point  radiators  on  the  platform  using  the 
CLEAN  algorithm.  Once  such  a  sparse  set  of  radiation  centers  have 
been  extracted,  it  is  possible  to  reconstruct  the  3-D  ACSAR  image  and 
the  3-D  frequency-spatial  data  very  rapidly  with  good  fidelity. 
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1  -  Introduction 

It  is  well  known  that  antenna  characteristics  can  be  dramatically  affected  by  the  presence 
of  platforms  that  support  the  antenna  structure.  Recently,  we  have  introduced  an  imaging 
algorithm  called  Antenna  Synthetic  Aperture  Radar  (ASAR)  that  images  the  locations  of  the 
dominant  scattering  off  the  platform  from  antenna  radiation  data  [1].  Such  information  is  useful 
to  the  antenna  designer  in  understanding  the  cause-and-effect  relationship  in  platform  scattering. 
In  this  work,  we  will  present  a  sparse  representation  of  platform  scattering  by  extracting  a  point 
scatterer  representation  of  the  ASAR  image.  Similar  to  an  ISAR  image,  an  ASAR  image  exhibits 
strong  point-scatterer  like  behavior.  In  the  ISAR  case,  these  points  are  called  scattering  centers 
[2].  In  our  case,  we  call  them  radiation  centers  since  they  are  in  fact  the  secondary  point 
radiators  on  the  platform.  First,  we  will  use  the  well-known  CLEAN  [3]  algorithm  to  extract  the 
radiation  centers.  Next,  we  will  show  that  we  can  reconstruct  and  extrapolate  the  antenna 
radiation  data  in  frequency  and  aspect  with  good  fidelity  by  using  the  radiation  centers.  Finally, 
we  will  demonstrate  that  the  radiation  center  information  can  be  used  to  provide  important 
guidelines  for  reducing  platform  effects. 

2  -  Extraction  of  Radiation  Centers 

Before  presenting  the  radiation  center  model  for  antenna  radiation  data,  we  shall  first 
review  the  ASAR  image  formation  algorithm  presented  earlier  in  [1]  and  provide  some 
motivations  for  the  radiation  center  model.  Using  the  Shooting  and  Bouncing  Ray  (SBR) 
technique  [4],  it  can  be  shown  that  a  3-D  ASAR  image  of  the  platform  can  be  generated  via  the 
expression: 

ASAR(x,y,z)  =  £Pj  htx-Xj.y-yj.z-Zj)  (1) 

i  rays 

where  h(x,  y,  z)  =  e--iko^r+x^  sjn  c(Ak  •  x)  -sin  c(k0A<)>  •  y)  •  sinc(k0  A9  •  z)  (2) 

In  the  above  expression,  (Xj,yj,Zj)  is  the  location  of  the  hit  point  on  the  target  for  the  ith  ray,  |5j 
is  the  ray  amplitude,  Ak,  A<|),  A0  are  the  half-bandwidths  in  the  k-,  <)>-  and  0 -domains, 
respectively.  A  key  observation  from  eq.  (1)  is  that  we  can  form  the  ASAR  image  using  SBR  on 
a  ray-by-ray  basis  directly  in  the  image  domain.  Each  ray  contributes  to  a  spot  centered  about 
(Xj,yj,Zj)  in  the  ASAR  image,  with  the  footprint  of  the  spot  governed  by  the  3-D  ray  spread 
function  h.  Since  tens  of  thousands  of  rays  are  traced  in  the  SBR  process,  this  may  imply  that 
tens  of  thousands  of  point  radiators  are  needed  to  adequately  represent  an  ASAR  image. 
However,  the  actual  ASAR  image  of  a  complex  platform  is  actually  quite  sparse.  This  is  because 
the  rays  interfere  with  one  another  to  give  rise  to  strong  coherent  scattering  over  only  a  small, 
localized  portion  of  the  target.  Therefore,  an  ASAR  image  can  be  accurately  represented  by  only 
a  limited  number  of  radiation  centers.  We  adopt  a  point  radiator  model  which  has  exactly  the 
same  form  as  (1)  to  parameterize  the  ASAR  data: 

ASAR(x,y,z)  =  XAnh(x-xn,y-yn,z-zn)  (3) 

n  rad  cntrs 

where  (x„,yn,zn)  is  the  radiation  center  location  and  An  is  its  amplitude.  To  determine  these 
parameters,  we  apply  the  image  processing  algorithm  CLEAN  which  has  already  been 
successfully  used  in  scattering  center  extraction  [2],  The  CLEAN  algorithm  iteratively  picks  out 
the  highest  point  in  the  image  and  removes  its  point-spread  response  from  the  image.  The 
extraction  process  is  iterated  until  the  maximum  in  the  residual  image  reaches  a  user-defined 


threshold.  As  an  example,  we  used  an  airplane  model  whose  CAD  file  (shown  in  Fig.l)  consists 
of  approximately  8,000  perfectly  conducting  facets.  A  dipole  antenna  is  located  above  the 
cockpit.  The  simulation  is  carried  out  with  a  SBR-based  code  Apatch  [5].  Fig.  3(a)  shows  the 
resulting  3-D  ASAR  images  at  different  z-slices.  The  center  frequency  is  10  GHz  with  0.56  GHz 
of  bandwidth  sampled  over  64  frequencies.  The  observation  angles  are  centered  about  the  nose 
of  the  airplane  and  range  between  -1.67°  and  1.67°  in  azimuth  sampled  over  32  points  and 
between  -0.86°  and  0.86°  in  elevation  over  8  points.  Fig.  3(b) -shows  the  locations  of  the 
extracted  radiation  centers  from  the  3-D  ASAR  image  of  Fig.  3(a)  using  the  3-D  CLEAN 
algorithm.  A  total  of  150  radiation  centers  are  extracted  from  the  image  and  their  locations  are 
plotted  in  Fig.  3(b)  as  small  circles.  We  notice  that  most  of  them  are  concentrated  around  the  z=- 
2m  cut  which  corresponds  to  the  top  surface  of  the  airplane.  The  rest  are  located  on  the  two  tails. 
Fig.  3(c)  shows  the  reconstructed  ASAR  image  from  the  150  extracted  radiation  centers.  The 
agreement  between  the  original  image  and  the  reconstructed  one  is  quite  good  over  a  dynamic 
range  of  50  dB. 

3  -  Applications: 

Two  examples  of  the  utilities  of  the  radiation  center  representation  are  presented.  In  the 
first  example,  we  demonstrate  how  the  radiation  center  concept  can  be  utilized  to  create  a  sparse 
model  of  the  antenna-platform  radiation  over  frequencies  and  aspect  angles,  thus  achieving  data 
compression.  To  fully  characterize  the  antenna-platform  radiation  pattern  for  all  angles,  we  need 
to  extract  3-D  radiation  center  models  at  various  angles  covering  the  entire  aspect  sphere.  In  our 
example,  we  focus  only  on  the  azimuth  angle.  The  platform-antenna  configuration  is  the  same  as 
before.  We  have  chosen  an  angular  granularity  of  5°  between  the  extraction  angles  of  each 
radiation  center  set  so  that  we  have  72  radiation  center  sets  in  total  covering  the  whole  360° 
azimuth.  Each  radiation  center  set  consists  of  100  radiation  centers.  The  bandwidths  used  for  the 
extraction  of  each  radiation  center  set  are  0.49  GHz  in  frequency  (with  a  10  GHz  center 
frequency),  2.22°  in  azimuth  and  2.46°  in  elevation.  We  have  reconstructed  the  azimuth  pattern 
from  0°  to  180°  by  using  72  radiation  center  sets  shown  as  the  dashed  line  in  Fig.  4(a).  To  test 
the  validity  of  the  model,  we  have  also  plotted  the  original  pattern  (drawn  as  solid  line) 
calculated  by  the  Apatch  code.  Since  the  antenna  pattern  on  a  complex  platform  is  such  a  rapidly 
varying  function  of  angle,  we  expand  the  30°-35°  region  in  Fig.  4(b)  to  better  display  the  fine 
agreement.  Note  thaF  the  data  between  30°-  32.5°  are  reconstructed  from  one  set  of  radiation 
centers  while  those  from  32.5°  -  35°  are  reconstructed  from  another  set.  The  data  between 
31.11°  and  33.89°  correspond  to  extrapolated  data  based  on  the  nearest-neighbor  radiation  center 
model.  The  overall  agreement  between  the  two  sets  of  data  are.,  fairly  good,  even  in  the 
extrapolated  region.  Therefore,  the  radiation  pattern  of  a  complex  antenna  pattern  can  be 
sparsely  represented  by  a  radiation  center  model  with  good  fidelity. 

In  the  second  example,  we  utilize  the  radiation  center  representation  to  carry  out  a 
platform  effect  reduction  study.  Since  the  radiation  centers  pinpoint  the  locations  of  the 
dominant  secondary  scattering  on  the  platform,  one  way  to  mitigate  such  undesirable  effects  is  to 
place  absorbers  at  those  locations.  Using  the  same  example  as  before,  we  have  identified  from 
the  database  those  radiation  centers  that  contribute  more  than  -25dB  to  the  radiated  field.  Then, 
by  tying  each  radiation  center  back  to  the  dominant  ray  hit  mechanism  that  gave  rise  to  it,  we 
place  403  absorbing  facets  around  the  platform.  The  absorber-coated  airplane  model  is  shown  in 
Fig.  2.  The  dark  triangles  on  the  platform  correspond  to  absorbing  facets.  Fig.  5  shows  the 
azimuth  radiation  pattern  from  the  platform  before  and  after  the  absorber  coating.  We  can  see 
that  for  most  of  the  angles,  the  scattered  field  from  the  platform  is  under  the  chosen  -25  dB  level. 
Therefore,  it  is  possible  to  suppress  platform  radiation  by  selectively  coating  the  platform  from 
the  radiation  center  information. 
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Figure  1.  View  of  the  platform  CAD  model  with  Figure  2.  The  CAD  model  of  the 

a  small  dipole  placed  above  cockpit.  absorber-coated  airplane. 


Figure  3.  (a)  3-D  ASAR  image  of  the  airplane  generated  from  the  simulated  data  about  nose 
using  the  Apatch  code,  (b)  The  location  of  the  1 50  radiation  centers  extracted  using 
CLEAN,  (c)  Reconstructed  3-D  ASAR  image  using  the  extracted  radiation  centers. 


(b) 


Figure  4.  Comparison  of  the  azimuth  radiation  patterns  from  the  original  data  (dashed)  and  the 
reconstructed  data  (solid)  based  on  radiation  centers,  (a)  Coarse  comparison  between 
0°  and  180°.  The  original  data  are  sampled  at  1°  while  the  reconstructed  one  are  sampled 
at  0.07°.  (b)  Detailed  comparison  between  30°  and  35°  at  the  same  angular  granularity. 


Figure  5.  Azimuth  radiation  patterns  before  (dashed)  and  after  (solid)  putting  the  absorbers  on  the 
platform. 
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1.  Introduction 

The  electromagnetic  coupling  between  antennas  on  a  complex  platform  is  an  important 
issue  in  antenna  design.  In  this  work,  we  set  out  to  develop  an  imaging  algorithm  to  pinpoint 
those  locations  on  the  platform  that  give  rise  to  antenna  interactions  based  on  the  coupling  data 
between  antennas.  Our  approach  to  this  problem  is  based  on  the  Inverse  Synthetic  Aperture 
Radar  (ISAR)  concept.  ISAR  imaging  is  a  standard  technique  to  map  the  locations  of  dominant 
scattering  off  a  target  based  on  the  multi-frequency,  multi-aspect  backscattered  data  [1],  We 
have  already  extended  this  concept  to  the  far-field  antenna  radiation  problem  by  introducing  the 
Antenna  SAR  (ASAR)  imaging  concept  [2].  By  collecting  the  multi-frequency,  multi-aspect 
radiation  data  from  an  antenna  in  the  presence  of  the  platform,  we  showed  that  it  is  possible  to 
image  the  locations  of  the  dominant  secondary  radiation  off  the  platform.  In  this  paper,  we  will 
further  extend  this  concept  to  the  near-field  antenna  coupling  problem  by  constructing  the 
Antenna  Coupling  SAR  (ACSAR)  image  of  the  platform.  We  will  show  that  a  platform  image 
can  be  generated  by  processing  multi-frequency,  multi-spatial  radiation  data  about  a  receiver 
location  from  a  transmitting  antenna.  Examples  of  the  ACSAR  image  will  be  presented  by  using 
the  computed  data  from  the  Shooting  and  Bouncing  Ray  (SBR)  code  Apatch  [3,4].  In  addition, 
we  will  present  a  fast  algorithm  for  ACSAR  image  formation  that  is  specifically  tailored  to  the 
SBR  technique.  By  taking  advantage  of  the  ray  tracing  information  within  the  SBR  engine,  we 
will  show  that  the  ACSAR  image  can  be  formed  with  little  additional  computation  time  than 
simple  ray  tracing. 

2.  ACSAR  Imaging  Algorithm 

We  assume  the  transmitter  and  receiver  setup  as  depicted  in  Fig.  1.  The  electric  field  at 
the  receiver  site  due  to  the  scattering  from  the  antenna-platform  interaction  at  point  P(xi,yi,zi) 
can  be  written  as 

\ 

Es(k)=  Aj  ■e"jkR"  e-jkR2i  (1) 

where  A.  is  the  strength  of  the  scattered  field,  Rn  is  the  path  length  from  transmitter  antenna  to  P, 
Rj,  is  the  path  length  from  P  to  the  receiver  location  and  k  is  the  free-space  wave  number.  Next, 
we  assume  that  the  radiation  data  are  collected  within  a  certain  frequency  bandwidth  and  on  a 
two-dimensional  aperture  around  the  receiver  point  (x0,  0,  0).  If  the  frequency  bandwidth  is 
small  compared  to  the  center  frequency  and  the  aperture  around  point  is  small  compared  to  the 
path  length  R^,  we  can  approximate  the  received  field  as  a  function  of  frequency  and  space  as 
follows: 


Es(k,y,z)  =  A;  •  e~jk(R'i+R2i)  .e-jk»TC“«'  ,e-jk0  *  “”P.  (2) 

In  the  above  formula,  it  is  obvious  that  there  is  a  Fourier  transform  relationship  between  the 
variables  (k,  y,  z)  and  (R=  Ru  +R2i  ,  u=k0cosa.  ,  v=k0sina,  sinp.  ).  Hence  by  taking  the  3-D 
inverse  Fourier  transform  of  the  scattered  field  with  respect  to  lc,  y,  and  z,  we  obtain  an  3-D 
image: 


ACSAR(R,u,v)=  Aj  •  8(R-Ri)-5(u-ui)-8(v- Vj) 


(3) 


Note  that  the  scattering  point  P  on  the  platform  will  manifest  itself  as  a  peak  in  the  ACSAR 
image  at  (R. ,  u. ,  v.)  with  amplitude  A..  It  is  straightforward  to  transform  the  image  from  the  (R, 
u,  v)  coordinates  to  the  desired  (x,  y,  z)  coordinates  of  the  platform.  Therefore,  by  taking  the  3- 
D  inverse  Fourier  transform  of  the  multi-frequency,  multi-spatial  scattered  field  data  about  the 
receiver  location  and  carrying  out  the  coordinate  transformation,  it  is  possible  to  construct  a  3-D 
image  of  the  platform  that  displays  the  strength  of  the  platform  scattering. 

As  an  example,  we  use  a  test  object  ‘slicy’  whose  CAD  geometry  is  shown  in  Fig.  2. 
The  platform  is  assumed  to  be  perfectly  conducting.  A  half-wave  dipole  is  used  as  the 
transmitter.  A  total  of  32  frequencies  are  computed  within  a  0.19  GHz  bandwidth  at  the  center 
frequency  of  10  GHz.  The  simulation  is  carried  out  by  an  SBR-based  code  Apatch  [4],  The 
scattered  field  is  collected  on  an  aperture  centered  at  (35,0,0).  There  are  32x32=1024  spatial 
points  on  the  aperture,  ranging  from  -0.49m  to  0.49m  in  z  and  from  -1.14m  to  1 .14m  in  y.  In  the 
computed  data,  we  include  only  the  scattered  field  from  the  platform  and  not  the  primary 
radiation  due  to  the  antenna  in  the  absence  of  the  platform.  The  resulting  3-D  ACSAR  image  is 
generated.  It  is  then  projected  into  three  2-D  projected  views  along  the  3  principal  planes  shown 
in  Fig.  3a.  We  observe  in  the  image  the  dominant  scattering  point  is  around  (17.5,  -40,  0)  which 
is  the  specular  point  between  the  transmitter  and  the  receiver.  We  also  observe  other  scattering 
mechanisms  on  the  top  surface  of  the  platform,  around  the  step  region  of  the  platform  and  at  the 
top  of  the  tall  cylinder. 


3.  Fast  ACSAR  Image  Formation  Using  SBR 

The  ACSAR  image  formation  outlined  in  the  last  section  requires  input  data  set  over 
multiple  frequencies  and  spatial  locations.  In  this  section,  we  will  derive  a  fast  imaging 
algorithm  specially  tailored  for  the  SBR  approach.  It  utilizes  the  ray  information  available 
within  the  ray  tracing  engine  to  generate  the  ACSAR  image  directly  on  a  ray-by-by  basis  without 
resorting  to  the  expensive  multi-frequency,  multi-spatial  calculations.  This  approach  is  similar 
to  the  fast  image  formation  algorithm  we  have  developed  previously  for  ISAR  and  ASAR  image 
simulation  [5,  6].  In  applying  the  SBR  technique  to  the  antenna  radiation  problem,  rays  are  first 
shot  from  the  phase  Renter  of  the  transmitter  and  traced  according  to  geometrical  optics.  At  the 
exit  point  of  each  ray  before  they  leave  the  platform  altogether,  a  ray-tube  integration  is  carried 
out  to  find  the  contribution  of  each  ray  to  the  near-field  at  the  receiver.  Under  this  construct,  the 
scattered  near  field  can  be  written  as; 

\ 

Es(k,  y,z)  =  •  e‘jk(R"+R2i>  •e-jk°ycos“i  •  e-jk»Z5ina>sinp‘  (4) 

i  rays 

This  equation  obtained  by  assuming  the  bandwidth  and  the  aperture  is  small  and  the  platform  is 
perfectly  conducting,  a,  is  proportional  to  the  field  at  the  exit  ray-tube,  and  is  only  weakly 
dependent  on  the  frequency  and  variations  on  the  path  length  R2i .  Substituting  (4)  in  to  our 
ACSAR  formula  and  performing  the  3D  inverse  Fourier  transform  operation,  we  get: 

ACSAR(R,u,v)=  IFT3{2ai  -e"jk(Rli+R2i)  •e'jlc<,ycosot|  .e-jk°*'sina‘sinp‘} 

i  rays 

=  £p;  hCR-Rj.u-Ui.v- Vj)  (5a) 

i  rays 

where 

h  (R,  u,  v)  =  e~jk°R  sin  c(  Ak  •  R )  •  sin  c(  Ay  •  u)  •  sin  c(  Az  •  v)  (5b) 

In  (5b),  Ak,  Ay,  Az  are  respectively  the  half  bandwidth  in  frequency  and  the  half  widths  of  the 
aperture.  We  can  rewrite  the  last  equation  as  a  convolution  between  a  3-D  weighted  impulse 
train  and  the  ray  spread  function  h: 

ACSAR(R,u,v)  =  (£pi-5(R-Ri,u-ui,v-vj))*h(R,u,v)=g*h 

i  rays 


(6) 


The  above  convolution  can  be  computed  quickly  using  an  interpolation  algorithm  and  FFT  [5,  6]. 

We  have  applied  the  fast  approach  to  the  same  geometry.  An  oversampling  ratio  of  4  is 
used  in  the  interpolation  grid  to  ensure  the  accuracy  in  the  computation.  The  results  are  depicted 
in  Fig.  3b  as  three  2-D  projected  ACSAR  images.  By  comparing  these  images  with  the  images 
that  are  generated  by  using  the  standard  frequency-aperture  approach,  we  observe  a  fairly  good 
agreement.  To  show  the  time  savings  that  we  have  gained  by  using  the  fast  algorithm,  we  list  the 
computation  times  in  Table  1.  The  ray-tracing  time  of  9  minutes  is  the  same  for  both  approach; 
however,  the  ray  summation  time  takes  about  106  hours  in  the  frequency-aperture  approach.  On 
the  other  hand,  the  fast  approach  reduces  the  ray  summation  time  is  to  only  6  minutes. 
Consequently,  the  combination  of  SBR  and  the  fast  ACSAR  technique  can  be  used  as  a  CAD 
tool  for  diagnosing  platform  effects. 
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Figure  1.  The  path  of  a  radiated  signal  that  is  scattered  off  a  point  on  the  platform. 
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I.  Introduction 

Moment  method  is  one  of  the  most  popular  frequency  domain  simulation 
methods  in  computational  electromagnetics  (CEM).  Since  it  solves  an  integral 
equation  rigorously  through  boundary  discretization,  huge  computer  memory  and 
computation  time  are  needed  to  simulate  the  scattering  phenomena  for  targets  of 
large  electrical  sizes.  Recently  the  development  of  the  fast  multipole  method 
(FMM)  and  the  multilevel  fast  multipole  algorithm  (MLFMA)  has  led  to  the 
reduction  of  the  computational  complexity  of  the  matrix-vector  multiplication 
operation  in  moment  method  from  DIN2)  to  <7(N,S)  or  D(NlogN)  [1,2].  The 
FLFMA-based  code  FISC  [3]  is  now  capable  of  computing  the  radar  cross  section 
(RCS)  of  a  full-size  fighter  airplane  at  1  GHz  on  a  workstation  in  several  hours  of 
run  time.  While  such  an  advance  in  computational  capability  has  brought  much 
excitement  to  the  CEM  community,  real-world  radar  signature  prediction  often 
requires  that  the  same  calculation  be  repeated  over  many  frequencies  and  aspect 
angles.  For  example,  in  automatic  target  identification  applications,  range  profiles 
of  the  target  are  commonly  used  to  pinpoint  target  features.  To  generate  range 
profiles  requires  hundreds  of  frequency  calculations  that  make  the  computation 
time  still  extraordinarily  long  even  using  FISC.  In  this  x paper,  we  present  a 
frequency  extrapolation  scheme  to  speed  up  the  signature  prediction  procedure 
using  FISC  by  avoiding  such  exhaustive  computations. 

We  adopt  a  model-based  approach  [4-6]  to  the  frequency  extrapolation 
problem.  The  key  in  the  success  of  an  extrapolation  algorithm  is  a  good  model  of 
the  physical  observable  to  be  extrapolated.  Here,  we  propose  a  scheme  that 
parameterizes  the  current  on  the  target,  which  is  available  in  the  FISC  solution, 
based  on  a  multipath  excitation  model.  We  first  run  FISC  for  the  scattering 
problem  at  several  frequency  points,  then  apply  the  ESPRIT  superresolution 
algorithm  [7]  directly  to  the  induce  current  output  at  these  frequency  points  to 
extract  the  time-of-arrival  and  amplitude  parameters.  Once  these  parameters  are 
obtained,  the  frequency  dependent  model  of  the  induce  current  on  each  facet  of 
the  target  can  be  constructed.  Thus  the  induced  currents  over  the  whole  frequency 
band  can  be  extrapolated  and  the  multifrequency  far  field  can  be  calculated.  The 
range  profiles  of  a  low-Q  missile  model  are  calculated  as  an  example  to 
demonstrate  the  performance  of  the  algorithm. 

II.  Model-Based  Parameterization 

We  assume  the  following  frequency  dependent  model  for  the  induced  current 
on  the  target: 


(1) 


J(f,S)  =  XAi(S)exp[-j27fti(S)] 

i=l 

where /is  the  frequency,  S  represents  the  position  on  the  target  surface,  and  (A,,  f.) 
are  respectively  the  amplitude  and  time-of-arrival  of  the  ith  scattering  mechanism 
observed  at  S.  This  is  a  physical  model  consistent*  with  high-frequency 
electromagnetic  phenomenology.  For  example,  for  a  point  on  the  target  surface  in 
the  lit  region  of  the  incident  wave,  the  First  term  of  the  expansion  should 
correspond  to  the  physical  optics  current.  The  higher  order  terms  correspond  to 
multiple  scattering  effects  such  as  currents  induced  by  multiple  reflected  or 
diffracted  fields  that  have  longer  time-of-arrival.  In  our  model,  we  assume  A.  to 
be  independent  of  frequency. 

With  this  model  in  hand,  we  next  extract  the  parameters  (A,,  t)  from  a  limited 
number  of  frequency  calculations  based  on  FISC  at  each  point  on  the  target.  We 
apply  the  superresolution  algorithm  ESPRIT  for  this  purpose  due  to  its  robustness 
in  the  presence  of  noise.  ESPRIT  is  based  on  the  data  model: 

F<v=S . /„  ra 

i~  1 

where  n(»)  denotes  additive  white  Gaussian  noise.  If  the  data  sequence  obeys  this 
ideal  model  and  the  number  of  the  sampling  points  M>2N+l,  ESPRIT  can 
estimate  N  and  resolve  each  A,  and  t, .  Of  course,  for  real  CEM  data,  the  nature  of 
the  noise  is  more  complex  and  will  degrade  the  performance  of  the  algorithm.  At 
frequencies  sufficient  above  target  resonances,  the  first  few  interactions  dominate 
the  scattering  phenomenon  and  only  a  small  N  is  needed  to  adequately  model  the 
observed  data.  Therefore  by  running  ESPRIT  using  the  FISC  generated  current  at 
M  frequencies,  the  frequency-dependent  current  model  in  (1)  can  be  obtained. 
Note  that  the  ESPRIT  processing  is  carried  out  at  each  point  on  the  target.  Once 
such  a  model  is  found,  we  can  extrapolate  the  data  to  determine  the  currents  at 
other  frequencies.  The  total  scattered  field  versus  frequency  can  thus  be 
calculated  by  integrating  the  extrapolated  currents. 

Two  comments  are  in  order.  First,  just  as  in  Fourier  processing,  the  maximum 
unambiguous  time  delays  that  can  be  resolved  by  ESPRIT  is  inversely 
proportional  to  the  frequency  sampling  of  the  data  points.  Therefore,  the 
frequency  sampling  must  be  chosen  sufficiently  small  to  avoid  aliasing  in  the 
time-of-arrival  estimates.  Typically,  we  choose  the  frequency  sampling  to  allow  a 
range  window  that  is  at  least  twice  the  size  of  the  target.  Second,  we  have 
adopted  the  time-of-arrival  model  for  the  induced  current  in  our  extrapolation 
scheme.  The  same  model  can  also  be  applied  directly  to  the  total  scattered  field. 
This  is  commonly  known  as  the  scattering  center  model.  However,  we  have 
found  that  the  scattered  field  requires  far  more  terms  to  model  adequately  for  a 
complex  target.  This  means  that  a  far  larger  number  of  frequency  points  are 
needed  from  FISC  simulation  to  achieve  a  good  frequency  extrapolation  model. 

III.  Numerical  Results 

To  demonstrate  the  performance  of  our  extrapolation  algorithm,  the  frequency 
response  and  range  profile  of  a  missile  model  are  computed  as  an  example.  The 
length  of  the  target  is  about  50  inches  and  the  model  consists  of  5679  perfectly 
conducting  facets  (see  Fig.  1).  The  brute-force  frequency  data  are  first  generated 
at  30°  in  azimuth  and  0°  in  elevation  over  a  frequency  band  from  2  GHz  to  6  GHz. 
The  frequency  sampling  interval  is  0.05  GHz  and  thus  a  total  of  81  frequency 


points  are  computed.  To  test  the  extrapolation  algorithm,  10  frequency  points 
between  2.50  GHz  and  2.95  GHz  from  FISC  are  used  to  extrapolate  the  scattered 
fields  to  other  frequencies  using  a  model  with  N=4.  The  frequency  response  of  the 
scattered  fields  for  the  VV  polarization  is  plotted  in  Fig.  2.  We  see  that  the 
extrapolated  far  fields  agree  quite  well  with  the  data  directly  computed  using 
FISC  within  the  frequency  range  2.5  to  4.5  GHz.  Discrepancies  in  amplitude 
occur  outside  of  this  range.  However,  note  that  the  peak  and  null  positions  of  the 
original  data  are  still  fairly  well  reproduced  by  the  extrapolation.  This  indicates  a 
good  time-of-arrival  estimate  and  the  discrepancy  is  more  likely  due  to  the 
frequency  dependent  amplitude  effects  which  are  not  included  in  the  present 
model.  Since  only  ten  frequency  points  near  the  low  frequency  end  are  computed, 
such  an  algorithm  leads  to  large  computational  savings.  The  total  computation 
time  (including  post-processing)  for  FISC  with  extrapolation  is  about  3  hours  in 
an  SGI  02  workstation  (R10K/155MHz),  while  it  takes  about  80  hours  on  the 
same  machine  for  the  brute-force  frequency  calculations.  Fig.  3  is  the  range 
profile  comparison  obtained  by  Fourier  transforming  the  frequency  response.  We 
see  that  the  main  scattering  features  of  the  extrapolated  results  coincide  very  well 
with  the  reference  results.  Fig.  4  shows  the  comparison  of  the  extrapolated  range 
profile  against  the  reference  result  for  the  HH  polarization.  Again,  good 
agreement  for  the  dominant  scattering  features  is  observed.  The  theoretical  and 
practical  limitation  of  this  extrapolation  algorithm  is  currently  under  study.  We 
will  also  report  on  the  extension  of  this  concept  in  simulating  two-dimensional 
radar  imagery  of  complex  targets. 
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Fig.l  The  outlook  of  the  missile  model. 


(inches) 

Fig. 3.  Comparison  between  the  range  profiles  for 
the  VV  polarization  obtained  by  Fourier 
transforming  the  frequency  response  in 
Fig.2. 
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Fig.2.  Comparison  of  the  backscattered  field 
vesus  frequency  between  the  brute-force 
FISC  results  and  the  extrapolation  results 
for  the  VV  polarization. 
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Fig.4.  Comparison  between  the  range  profiles  for 
the  HH  polarization  obtained  by  brute- 
force  FISC  and  FISC  +  extrapolation. 
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I.  Introduction 

The  impedance  boundary  condition  (IBC)  approximation  is  an  effective  way 
to  model  complex  structures  [1-3].  It  replaces  the  original  structure  with  a  surface 
impedance  so  that  the  problem  dimension  is  reduced  by  one.  In  large-scale 
problems  such  as  the  scattering  from  a  3-D  target  coated  with  absorbing  materials, 
the  IBC  approach  will  result  in  large  savings  in  computer  resources  since  the  fine 
details  of  the  structure  are  replaced  by  simple  impedance  boundary  conditions  on 
the  surface.  In  this  paper,  we  set  out  to  determine  the  equivalent  IBC  model  of  a 
material-coated  corrugated  grating  (Fig.  1). 

It  is  not  an  easy  task  to  find  a  simple  EBC  that  accurately  approximates  the 
original  structure  over  a  broad  range  of  incident  angles  for  an  arbitrary  grating 
structure.  We  therefore  employ  the  genetic  algorithm  (GA)  to  find  the  optimal 
parameters  in  our  EBC  model.  In  our  approach,  a  rigorous  solution  of  the 
reflection  coefficients  at  a  few  incident  angles  is  first  calculated  using  a  method  of 
moments  (MoM)  code.  The  IBC  model  is  used  to  predict  the  reflection 
coefficients  at  the  same  observation  angles.  The  model  coefficients  are  then 
optimized  so  that  the  difference  between  the  IBC-predicted  and  the  MoM- 
predicted  reflection  coefficients  is  minimized.  The  genetic  algorithm  proves  to  be 
quite  efficient  in  obtaining  an  accurate  IBC  model.  The  agreement  between  the 
IBC-predicted  reflection  coefficients  and  the  reference  MoM  ones  are  good  at 
most  incident  angles,  even  near  grazing. 

II.  Formulation 

In  our  approach,  the  reflection  coefficients  at  a  number  of  observation  angles 
are  first  computed  by  the  method  of  moments.  This  solution  serves  as  the 
reference  data  for  our  EBC  model.  The  MoM  formulation  for  a  material  coated 
grating,  as  shown  in  Fig.  1,  begins  by  dividing  one  cell  of  the  grating  into 
different  homogenous  regions  according  to  the  material  layers.  The  boundary 
integral  equations  are  then  obtained  for  each  region.  Periodic  boundary  conditions 
are  applied  in  each  region  and  the  fields  are  matched  at  the  interfaces  between 
regions.  Finally,  the  fields  in  the  top  half-space  are  expanded  into  a  sum  of 
Floquet  harmonics  and  are  matched  to  the  fields  in  the  lower  region  so  that  the 
reflection  coefficients  can  be  calculated.  A  detailed  description  of  this  method 
can  be  found  in  [4]. 

Next,  we  determine  the  reflection  coefficients  based  on  an  approximate  IBC 
model.  The  impedance  boundary  condition  can  be  written  as  [1] 
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Since  the  structure  is  periodic,  each  element  of  the  boundary  impedance  can  be 
expanded  in  a  Fourier  series  and  the  expansion  coefficients  are  the  model 
parameters  to  be  determined.  Note  that  the  cross  impedance  terms  Zu  and  Za  are 
kept  in  our  formulation  so  we  can  assess  their  importance.  The  reflected  wave 
from  the  periodic  structure  consists  of  an  infinite  number  of  Floquet  harmonics, 
either  propagating  or  evanescent.  By  truncating  the  Floquet  harmonics  to  a  finite 
order  and  matching  fields  at  the  boundary,  a  matrix  equation  can  be  set  up  for 
determining  the  reflection  coefficients.  This  relationship  is  utilized  by  the  genetic 
algorithm  to  evaluate  the  IBC  parameters. 

In  the  final  step,  the  model  parameters  are  optimized  by  the  genetic  algorithm. 
In  GA,  the  parameters  to  be  optimized  are  first  encoded  into  binary  form  [5,  6].  A 
set  of  the  encoded  parameters  is  known  as  a  chromosome.  The  basic  idea  of  GA  is 
to  generate  a  pool  of  chromosomes,  discard  the  bad  ones,  keep  the  good  ones  and 
let  them  evolve  to  produce  better  chromosomes.  Each  searched  chromosome  is 
evaluated  by  a  cost  function  which,  in  our  case,  is  the  mean  square  error  between 
the  reflection  coefficients  computed  by  MoM  and  those  predicted  by  the  IBC 
model.  The  evolution  process  is  repeated  until  the  best  parameters  are  found. 

Compared  with  other  IBC  approaches,  the  method  described  above  has  several 
advantages.  First,  the  boundary  impedance  is  assumed  to  be  anisotropic  so  that  the 
same  model  can  be  applied  to  oblique  incidence  from  any  angle  0  and  <j>  shown  in 
Fig.  1.  Second,  it  is  possible  to  build  in  spatial  variation  of  the  boundary 
impedance  by  adjusting  the  number  of  terms  used  in  the  Fourier  series  expansion. 
By  using  moreTerms,  the  IBC  model  can  be  made  more  accurate.  Moreover,  the 
boundary  impedance  is  solved  from  error  minimization  by  the  genetic  algorithm. 
Thus  the  number  of  observation  points  is  flexible  and  can  be  increased  to  obtain 
better  results.  Most  important  of  all,  the  position  of  the  impedance  surface  can 
also  be  optimized.  In  general  layered  coating  problems,  there  is  not  a  preferred 
position  for  the  impedance  surface.  The  solution  can  sometimes  be  improved  by 
applying  the  IBC  at  a  position  other  than  a  natural  interface  in  the  structure  [3]. 
Therefore,  by  solving  for  the  best  position  of  the  impedance  surface  as  one  of  the 
model  parameters,  the  accuracy  of  the  model  can  be  improved. 

III.  Results 

In  this  section,  two  examples  are  presented  to  demonstrate  the  effectiveness  of 
the  method.  The  first  example  is  a  deep  triangular  grooved  grating  with  relatively 
small  period.  The  geometry  of  one  cell  of  the  grating  is  shown  in  Fig.  2  (a).  The 
coating  material  is  MagRAM  with  the  material  constants  e  =  14.35  - y'0.28  and  nr 
=  1.525  -  y  1.347  at  the  frequency  of  10  GHz.  17  observation  angles  are  selected 
which  include  normal  incidence  and  the  combination  of  0  =  20°,  40°,  60°,  80°  and 
<p  =  0°,  30°,  60°,  90°.  In  order  to  compare  with  the  simple  tensor  impedance 
boundary  condition  (TIBC)  approach,  we  set  Za  and  in  (1)  to  zero  and  use  only 
one  term  in  the  Fourier  expansion  for  Za  and  Za.  The  reflection  coefficients  of  the 
H-pol  and  V-pol  incidence  are  plotted  in  Figs.  3  (a)  and  3  (b),  respectively.  In  the 
figures,  the  x  axis  is  divided  into  sections  of  different  incident  angle  0,  which 
varies  from  5  to  85  degrees  in  steps  of  10  degrees.  In  each  of  the  0  section,  the 


grating  angle  <p  varies  from  0  to  90  degrees.  In  using  the  TEBC,  we  put  the 
impedance  surface  at  the  conductor  backing  since  it  is  the  only  plane  in  the 
structure.  It  is  obvious  that  TIBC  is  good  only  when  the  incident  angle  9  is  small, 
especially  for  the  V-pol  incidence.  But  with  the  GA  approach,  the  position  of  the 
impedance  surface  is  also  optimized  so  that  the  matching  is  good  even  near 
grazing  incidence. 

In  example  2,  we  compare  the  IBC  approximations  with  and  without  the  cross 
impedance  terms  and  Za.  The  geometry  of  the  structure  is  shown  in  Fig.  2  (b). 
The  material  constants  are  e  =  8.3-y'2.4  and  pr  =  2-j0.9  at  10  GHz.  The  same 
observation  points  are  used  as  in  the  previous  example.  The  0th  order  IBC  is 
determined  and  the  comparison  between  the  reflection  coefficients  is  illustrated  in 
Fig.  4.  With  the  cross  impedance  included,  the  accuracy  of  the  approximation  is 
improved.  Further  improvement  in  the  accuracy  can  be  achieved  by  increasing 
the  order  of  the  IBC  model. 

IV.  Conclusion 

In  this  paper,  an  impedance  boundary  condition  model  based  on  the  genetic 
algorithm  is  proposed  to  approximate  arbitrary  coated  gratings  in  scattering 
problems.  The  periodic  structure  is  replaced  by  a  periodic  IBC  on  a  virtual 
surface.  The  boundary  impedance  and  the  position  of  the  surface  are  optimized  by 
matching  the  reflection  coefficients  to  the  rigorous  numerical  solution  at  a  number 
of  incident  angles.  Similar  to  traditional  EBC  models,  this  approach  is  most 
effective  when  the  coating  material  is  high-loss  and  of  high  contrast.  Finally,  the 
output  IBC  model  generated  by  this  algorithm  can  be  incorporated  into  an  existing 
computational  electromagnetics  code  to  assess  the  performance  of  the  corrugated 
surface  when  mounted  on  complex  platforms. 
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I.  Introduction 

In  antenna  design  and  analysis,  the  mounting  platform  can  have  a  significant  effect  on  the 
antenna  radiation  characteristics.  However,  rigorous  solution  of  the  radiation  problem  over  a  complex 
platform  is  very  time  consuming,  and  the  computation  complexity  increases  dramatically  as  the 
frequency  increases.  In  this  paper,  we  present  a  model-based  frequency  extrapolation  technique  with 
which  the  radiated  field  over  a  broad  band  of  frequencies  can  be  obtained  using  the  rigorously 
computed  results  at  low  frequencies. 

Our  approach  entails  three  steps.  First,  the  induced  current  on  the  platform  surface  is  computed 
at  low  frequencies  using  the  method  of  moments  (MoM).  Second,  we  apply  the  time-of-arrival  model 
to  the  current  on  each  basis  element  on  the  surface.  The  model  coefficients  are  obtained  using  super¬ 
resolution  algorithm  ESPRIT  [1].  Finally,  the  induced  current  at  higher  frequencies  is  computed  using 
the  model  and  the  radiation  characteristics  are  calculated.  This  approach  is  similar  to  that  present  in  [2] 
and  [3]  for  radar  signature  extrapolation.  Our  results  show  that  when  the  frequencies  and  the 
discretization  of  the  platform  are  properly  chosen,  the  radiated  field  at  higher  frequencies  can  be 
extrapolated  with  only  moderate  computational  cost. 

\ 

II.  Extrapolation  Methodology 

As  the  first  step  of  the  extrapolation  process,  the  induced  current  on  the  target  is  computed  at  a 
small  set  of  points  at  low  frequencies.  Once  the  current  over  the  surface  is  computed  at  low 
frequencies,  we  apply  the  time-of-arrival  model  to  each  of  the  current  element.  In  this  model,  we 
assume  that  the  total  current  is  induced  by  different  scattering  mechanisms,  as  shown  in  Fig.  1 .  Each 
of  the  incident  mechanism  has  a  distinct  arrival  time,  so  that  the  current  can  be  written  as 

J(0))  =  fjane-j^  (1) 

n=l 

where  co  is  the  angular  frequency,  N  is  the  total  number  of  incident  mechanisms  and  t„  is  the  arrival 
time  of  the  incident  mechanism. 

Since  the  incident  mechanisms  correspond  to  scattering  from  different  parts  of  the  target,  the 
maximum  difference  in  tn  is  related  to  the  size  of  the  target.  Thus  the  sampling  rate  in  the  frequency 
domain  should  be  high  enough  to  distinguish  these  time  events  from  all  parts  of  the  target.  Based  on 
this  consideration,  we  approximately  constraint  the  sampling  rate  in  frequency  to 


(2) 


Af<c/D 

where  c  is  the  speed  of  light  and  D  is  the  maximum  dimension  of  the  target. 

The  model  coefficients  an  and  tn  are  obtained  using  the  superresolution  algorithm  ESPRIT, 
which  is  based  on  the  model  that  the  signal  consists  of  a  sum  of  exponential  and  additive  white  noise. 
Given  a  sequence  with  M  samples,  the  algorithm  can  estimate  the  number  of  exponential  N  and 
determines  the  amplitude  and  period  of  each  exponential  term.  The  basic  requirement  in  the  number  of 
samples  is  M  >  2N  +  1.  Once  the  model  parameters  are  found,  the  induced  current  at  higher 
frequencies  can  be  computed  using  (1).  The  radiated  field  is  then  easily  obtained  from  the  extrapolated 
current. 


III.  Results 

As  an  example,  we  consider  a  2-D  structure  as  shown  in  Fig.  2  (a).  We  are  interested  in  the 
radiation  pattern  over  a  frequency  band  from  0.15  to  0.45  GHz.  To  obtain  the  data  for  the 
extrapolation,  we  compute  the  induced  current  at  10  frequencies  from  0.15  to  0.24  GHz.  Then  we  use 
the  ESPRIT  algorithm  to  obtain  the  model  coefficients  for  each  current  element  and  compute  the 
radiated  field  based  on  the  model.  Fig.  2(a)  shows  the  radiated  field  of  a  horizontally  polarized  line 
source  as  a  function  of  frequency  at  an  elevation  angle  of  40°.  The  dashed  curve  is  obtained  from  the 
model-based  extrapolation  while  the  solid  curve  is  the  reference  brute-force  solution.  We  observe  that 
the  extrapolation  algorithm  correctly  predicts  the  peaks  and  null  positions  in  frequency,  indicating  a 
good  estimate  on  the  times-of-arrival.  Fig.  2(b)  is  the  time  domain  response  obtained  via  an  inverse 
Fourier  transform  of  the  frequency  data.  The  first  large  peak  corresponds  to  the  specular  scattering 
from  the  flat  surface  and  the  second  large  peak  is  due  to  the  scattering  from  the  step  region. 

We  notice  from  Fig.  2(a)  that  the  model-predicted  field  matches  well  with  the  reference  field  at 
the  first  10  frequencies,  but  drops  below  the  computed  field  as  frequency  goes  higher.  This  is  because 
the  intensity  of  the  field  radiated  by  the  line  source  is  proportional  to  the  square  root  of  the  frequency. 
We  can  overcome  this  by  compensating  this  effect  before  doing  the  extrapolation.  Thus  the  time-of- 
arrival  model  becomes 


(3) 

vn=1  J 

where  4a  should  be  replaced  by  co for  3-D  situations.  After  the  compensation,  the  extrapolation  result 
is  further  improved,  as  shown  by  the  frequency  and  time  responses  in  Figs.  3(a)  and  3(b). 


Finally  we  look  at  the  radiation  problem  of  the  3-D  platform  in  Fig.  4(a).  The  source  is  a 
horizontally  polarized  dipole.  The  solver  used  is  FISC  [4],  which  is  a  3D  MoM  code  based  on  the  fast 
multipole  method.  Similar  to  the  2-D  case,  we  use  the  computed  current  at  10  frequencies  from  0.15  to 
0.24  GHz  to  extrapolate  the  data  to  0.45  GHz.  The  extrapolated  radiation  field  as  a  function  of 
frequency  at  the  elevation  angle  of  40°  is  shown  in  Fig.  4(b)  as  the  dashed  curve.  The  reference  brute- 
force  solution  is  plotted  as  the  solid  curve.  The  major  features  of  the  radiated  field  are  well 
characterized  by  the  extrapolation  algorithm. 


IV.  Conclusion  and  Discussion 


As  we  have  seen  from  the  results,  the  frequency  extrapolation  technique  is  an  efficient  way  of 
obtaining  the  radiation  pattern  over  a  broad  band  of  frequencies.  Computation  time  is  reduced 
dramatically  since  the  current  is  rigorously  solved  only  at  low  frequencies.  We  improved  the  result  by 
compensating  the  frequency  factor  of  the  source  in  the  time-of-arrival  model.  This  indicates  that  a 
wrong  frequency  dependence  in  the  model  may  result  in  errors  in  the  model  coefficients.  In  addition  to 
the  frequency  factor  of  the  source,  the  frequency  dependence  of  different  scattering  components  of 
each  current  element  could  be  different,  due  to  different  scattering  physics.  Further  incorporation  of 
these  effects  should  further  enhance  the  accuracy  of  the  extrapolation. 
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Fig.  3.  Frequency  compensated  extrapolation  (H-pol). 
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I.  Introduction 

The  characterization  of  radiation  from  a  complex  platform  is  usually  very  time 
consuming,  since  the  computation  complexity  of  the  rigorous  solution  of  a 
radiation  problem  grows  dramatically  as  frequency  goes  high.  On  the  other  hand, 
a  broad  bandwidth  and  a  high  center  frequency  are  usually  required  in  order  to 
characterize  the  radiation  features  of  the  platform.  In  this  paper,  we  present  a 
frequency  extrapolation  technique  for  the  antenna  radiation  problem  to  obtain  the 
radiated  field  over  a  wide  band  of  frequencies.  A  radiation  center  representation 
of  the  platform  is  then  found  from  the  extrapolated  data. 

We  start  with  a  CAD  model  of  the  platform.  First,  the  induced  current  on  the 
surface  is  computed  using  the  method  of  moments  (MoM)  at  a  set  of  low 
frequencies.  A  time-of-arrival  model  is  then  applied  to  each  current  element  and 
the  super-resolution  algorithm  ESPRIT  [1]  is  used  to  calculate  the  model 
coefficients.  The  frequency  response  of  the  current  is  extrapolated  based  on  this 
model  and  the  radiated  field  is  obtained  over  a  wide  band  of  frequencies.  Next,  we 
assume  the  radiation  can  be  approximated  by  the  radiation  of  a  set  of  radiation 
centers.  We  estimate  the  locations  of  the  radiation  centers  by  generating  a 
Antenna  Synthetic  Aperture  Radar  (ASAR)  [2]  image,  which  is  simply  an  Inverse 
Fourier  Transform  of  the  radiated  field.  Fine  positions  of  each  radiation  center  is 
found  using  match  pursuit  technique  and  the  radiation  centers  are  extracted  one 
after  another  until  the  remaining  energy  is  small  enough  to  be  ignored.  Thus  we 
come  to  a  radiation  center  representation  of  the  original  radiation  problem. 

Results  show  that  the  frequency  extrapolation  is  quite  accurate  and  the 
computation  time  is  saved  considerably.  The  radiation  center  model  characterizes 
most  radiation  features  of  the  platform.  It  gives  a  sparse  representation  of  the 
radiation  problem  and  is  very  useful  in  analysis  and  design  problems. 

II.  Methodology 

At  the  beginning  of  the  process,  we  compute  the  induced  current  at  a  few 
low  frequencies  using  MoM,  using  a  CAD  model  as  the  one  in  Fig.  1.  We  assume 
that  the  total  current  is  induced  by  different  scattering  mechanisms  [3]  [4].  Each  of 
the  incident  mechanism  has  a  distinct  arrival  time.  Also  taking  into  account  the 
fact  that  the  amplitude  of  the  radiated  field  from  a  dipole  source  is  proportional  to 
the  frequency,  a  current  element  can  be  written  as 


(1) 


J(co)  =  <o[]  ane~i<atn 

rtm  1 

where  co  is  the  angular  frequency,  N  is  the  total  number  of  incident  mechanisms, 
a„  and  t„  are  the  amplitude  and  arrival  time  of  the  /Ith  incident  mechanism, 
respectively. 

The  model  coefficients  a„  and  t„  are  obtained  using  ESPRIT,  which  is  based 
on  the  model  that  the  signal  consists  of  a  sum  of  exponential  and  additive  white 
noise.  Once  the  model  parameters  are  found,  the  induced  current  at  higher 
frequencies  can  be  computed  using  (1).  The  radiated  field  is  then  easily  obtained 
from  the  extrapolated  current. 


We  assume  that  the  radiated  field  can  be  approximated  by  the  radiation  of  a 
set  of  point  radiators.  To  estimate  the  location  of  the  radiation  centers,  we  Inverse 
Fourier  Transform  the  radiated  field  with  respect  to  the  frequency  and  incident 
angles  to  generate  an  ASAR  image.  The  detail  of  the  ASAR  algorithm  can  be 
found  in  [2],  The  point  with  the  highest  intensity  is  located  in  the  image  and  its 
amplitude  and  coordinate  are  determined  to  serve  as  an  estimation  of  the  strongest 
radiation  center.  Then  we  perform  a  match  pursuit  algorithm  to  find  the  exact 
location  [5].  Assume  a  radiation  center  is  close  to  the  estimated  point,  we  can  find 
the  radiated  field  from  the  radiation  center  given  its  location  and  the  position  of 
the  antenna  by 


E\k,dx  ,02 )  -  Ae~Jkr°e 


jk(xQ  sin  9  cos  j+y0  sin  6  sin^ +z0  sin  6  ) 


M  g>(***°+ Vo  +M. ) 
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where  k  is  the  wave  number,  (*&  yo,  zo)  is  the  location  of  the  radiation  center.  The 
angles  9i,  62  are  defined  in  Fig.  2.  By  slightly  varying  the  position  of  the  radiation 
center  and  subtracting  its  radiated  field  from  the  total  field,  we  can  come  up  with 
an  optimal  location  that  causes  the  maximum  energy  reduction.  Then,  subtract  the 
contribution  of  this  radiation  center  from  the  total  radiated  field  and  repeat  the 
above  process,  we  are  able  to  obtain  a  set  of  radiation  centers,  until  the  remaining 
energy  is  small  enough. 


III.  Results 

An  example  of  the  antenna-platform  problem  is  shown  in  Fig.  1.  The  induced 
current  is  rigorously  computed  at  the  frequencies  from  0.15  to  0.24  GHz  at  an 
interval  of  0.01  GHz.  The  computation  is  done  by  FISC  [6],  which  is  a  MoM  code 
based  on  the  fast  multipole  method  (FMM).  The  current  is  then  extrapolated  and 
the  radiated  field  is  computed.  The  extrapolation  result  is  compared  with  the 
reference  field  data,  which  is  also  computed  by  FISC  over  the  frequency  band 
from  0.15  to  0.45  GHz  at  6  -  50°  and  ^  =  0°,  as  shown  in  Fig.  3(a).  The 
corresponding  time  domain  response  is  plotted  in  Fig.  3  (b).  The  matching 
between  the  two  is  good  in  both  the  frequency  and  time  domain. 

Next,  we  extract  the  radiation  centers  using  the  match  pursuit  algorithm  base 
on  the  extrapolated  field  over  the  range  that/ varies  from  0.6  to  0.9GHz.  The  first 


10  extracted  radiation  centers  are  shown  in  Fig.  5.  As  we  can  observe,  most  of  the 
radiation  centers  are  located  on  the  surface  of  the  platform  as  we  expected. 

The  radiated  field  can  be  easily  reconstructed  using  the  radiation  centers.  The 
original  radiated  field  at  kz  =  n  is  plotted  in  Fig  6  (a)  as  a  function  of  the  kx  and  ky, 
where  kx,  ky  and  kz  are  defined  in  (2).  The  field  reconstructed  from  the  first  10 
radiation  centers  is  shown  in  Fig  6  (b).  The  two  patterns  match  well  at  most  of  the 
incident  angles.  Thus  the  radiation  pattern  is  efficiently  approximated  by  the 
radiation  center  model. 

IV.  Conclusions 

We  have  applied  the  frequency  extrapolation  technique  to  the  antenna 
radiation  problem.  The  radiated  field  over  a  wide  band  of  frequencies  can  be 
extrapolated  efficiently  based  on  the  computation  result  at  a  few  lower 
frequencies.  A  radiation  center  model  of  the  platform  can  be  found  by  performing 
the  match  pursuit  algorithm.  The  model  proves  effective  in  approximating  the 
original  problem  since  the  radiated  field  can  be  accurately  reconstructed,  with 
very  little  computation  effort. 
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Fig.  1.  Geometry  of  3-D  platform 


(a)  Frequency  domain  response  (dB) 


Fig.  3.  Frequency  Extrapolation 
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Fig.  2.  Radiation  center  model 
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Fig.  4.  Extracted  radiation  centers 
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(a)  Original  radiated  field 
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(b)  Reconstructed  radiated  field 
Fig.  6.  Reconstruction  of  the  radiated 
field  from  the  radiation  center  model 


